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The  techniques  evaluated  included  two  regional  methods,  Midwest  Research 
Institute  (MRI)  and  National  Eutrophication  Survey  (NES).  The  three  instream 
techniques  evaluated  were_^  average  flow,  average  concentration  (Q*C);  flow- 
weighted  concentration  (QC);  and  flow  interval  (FI).  Regional  methods  re¬ 
quired  delineation  of  watershed  land  use,  soil  descriptions,  and  consultation 
with  experts,  e.g.,  U.  S.  Department  of  Agriculture-Soil  Conservation  Service, 
for  coefficient  selection.  Instream  techniques  required  streamflow  and  water 
quality  data. 

The  agricultural  Sandusky  River  and  Honey  Creek  watersheds  in  Ohio  and 
the  forested  Caddo  River  watershed  in  Arkansas  were  selected  as  prototypes  for 
comparison  of  these  techniques.  Water  quality  parameters  included  in  the  study 
were  SS,  TKN,  and  TP.  The  results  were  used  as  input  to  a  reservoir  trophic 
status  assessment  procedure  to  compare  their  applicability  to  Corps  of  Engi¬ 
neers  (CE)  reservoir  planning. 

As  a  result  of  this  study,  it  was  concluded  that  the  primary  value  of  the 
MRI  technique  was  to  familiarize  the  user  with  the  watershed’s  physiographic 
features  and  hydrologic  response  characteristics  to  better  interpret  water 
quality  loading  estimates.  The  NES  method  has  only  limited  application,  at 
best,  to  CE  projects  in  selected  geographical  areas  of  the  United  States.  The 
generalized  regional  methods  for  the  most  part  were  not  adequate  for  estimat¬ 
ing  annual  water  quality  loadings  to  reservoirs. 

Of  the  three  techniques  that  use  available  surface  water  quality  and 
streamflow  data,  only  the  FI  technique  is  statistically  valid,  and  it  is  recom¬ 
mended  for  estimating  annual  water  quality  loadings  to  reservoirs  based  on  the 
aforementioned  evaluation  factors.  However,  those  water  quality  parameters 
that  may  load  streams  as  a  function  of  season  and  flow  would  require  intensive 
sampling  and  analyses  to_statistically  characterize  the  short-term  loadings 
during  the  ^ear.  The  Q*C  technique  is  applicable  for  water  quality  consti¬ 
tuents  that  maintain  fairly  constant  concentrations  through  the  range  of 
streamflow.  The  QC  technique  requires  data  representative  of  the  flow  range 
and  runoff  water  quality  characteristic  of  the  watershed. 

Field  sampling  programs  for  the  purpose  of  estimating  nonpoint  source 
water  quality  loadings  to  reservoirs  should  include  simultaneous  measurement  of 
flow  and  water  quality  to  the  reservoir.  Samples  should  represent  wide  ranges 
of  flow  at  different  seasons  of  the  year.  Automatic  water  samplers  that  are 
activated  by  changes  in  flow  velocity  or  river  stage  are  recommended  where 
extensive  sampling  is  required. 
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mental  Laboratory  (EL). 
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Dr.  John  Harrison,  Chief,  EL.  The  report  was  written  by  Dr.  Westerdahl, 
Mr.  Ford,  Ms.  Harris,  and  Dr.  Lee. 

Appreciation  is  expressed  for  manuscript  review  to  Mr.  L.  Mulkey, 
Environmental  Protection  Agency,  Athens,  Georgia;  Dr.  M.  Frere,  U.  S. 
Department  of  Agriculture  (USDA) -Southern  Regional  Office,  New  Orleans, 
Louisiana;  Dr.  D.  DeCoursey,  USDA-Sedimentation  Laboratory,  Oxford,  Mis¬ 
sissippi;  Mr.  E.  Lawson,  U.  S.  Forest  Service,  Fayetteville,  Arkansas; 
Mr.  S.  Ursic,  U.  S.  Forest  Service,  Oxford,  Mississippi;  Messrs.  M. 
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EVALUATION  OF  TECHNIQUES  TO  ESTIMATE  ANNUAL 
WATER  QUALITY  LOADINGS  TO  RESERVOIRS 

PART  I:  INTRODUCTION 
Background 

1.  Surface  runoff  results  in  nonpoint  source  loadings  of  chemical 
constituents  and  suspended  solids,  which  directly  influence  reservoir 
project  purposes,  management  operations,  and  reservoir  eutrophication. 
Over  an  extended  period  of  time,  the  chemical  and  sediment  contribution 
to  the  project  from  surface  runoff  and  baseflow  may  change  significantly 
as  a  result  of  the  gradual  process  of  watershed  urbanization  or  major 
land  use  modifications.  Hence,  short-  and  long-term  influence  of  water¬ 
shed  runoff  on  existing  and  proposed  projects  must  be  understood  and 
predictable . 

2.  With  existing  technology,  field  personnel  frequently  have  dif¬ 
ficulty  accurately  assessing  the  influence  of  nonpoint  source  chemical 
and  sediment  loadings  on  the  water  quality  of  existing  and  proposed  res¬ 
ervoir  projects.  Sufficient  data,  funding,  and  time  usually  are  not 
available  in  early  planning  stages  to  allow  Corps  of  Engineer  (CE)  per- 
lonnel  to  apply  existing  simulation  models,  e.g.  Hydrocomp  Interna¬ 
tional,  Inc.  (1976)  and  U.  S.  Army  Hydrologic  Engineering  Center  (1977). 

3.  The  requirement  for  predicting  water  quality  in  proposed  res¬ 
ervoirs,  having  been  established  by  previous  Public  Laws  and  Executive 
Orders,  has  contributed  to  the  development  of  several  simplified  tech¬ 
niques.  These  techniques  require  evaluation  to  determine  which  one 
provides  estimates  of  annual  loadings  to  reservoirs  suitable  for  predict¬ 
ing  long-term  effects  of  nonpoint  chemical  and  sediment  loadings  on  res¬ 
ervoir  water  quality  and  eutrophication. 

Purpose 

4.  The  purpose  of  this  study  was  to  evaluate  techniques  most 
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commonly  used  by  CE  District  personnel  to  estimate  mean  annual  export 
rates  of  nonpoint  source  water  quality  constituents  from  watersheds  of 
proposed  and  existing  reservoirs. 


Scope 

5.  This  report  presents  comparative  predictions  of  annual  water 
quality  loadings  to  reservoirs  using  techniques,  exclusive  of  hydrologic 
and  water  quality  modeling,  most  commonly  used  for  planning  by  the  CE 
District  Offices.  The  techniques  were  used  as  they  exist;  no  modifica¬ 
tions  or  improvements  were  made.  Since  comparisons  of  techniques  were 
made  using  existing  field  data,  prototype  watersheds  were  selected  for 
which  extensive  stream  water  quality  and  streamflow  data  were  readily 
available,  including  detailed  land  use,  soil  type,  and  resource  manage¬ 
ment  information.  These  included:  (a)  Caddo  River  watershed  above 
Glenwood,  Arkansas;  (b)  Sandusky  River  watershed  above  Fremont,  Ohio; 
and  (c)  Honey  Creek  watershed  above  Melmore,  Ohio.  Detailed  discussion 
of  field  sampling  strategies  was  not  within  the  scope  of  this  report. 
Emphasis  was  placed  on  the  application  of  each  technique  and  the  reason¬ 
ableness  of  annual  water  quality  loading  predictions  derived  from  the 
same  field  data  and  land  use  information  for  each  prototype  watershed. 

Approach 

6.  After  a  review  of  numerous  approaches  for  estimating  loadings, 
e.g.  Lystrom  et  al.  (1978),  True  (1976),  Wu  and  Ahlert  (1978),  Johnson 
et  al.  (1976),  Fisher  et  al.  (1968),  and  Biggar  and  Corey  (1969),  it  was 
concluded  that  an  interagency  workshop  should  be  held  to  select  the  ap¬ 
proaches  to  be  evaluated.  In  February  1978,  this  workshop  was  convened 
at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg, 
Mississippi.  Recommendations  from  the  participants  included  evaluation 
of  existing  techniques  and  data  requirements  to  predict  annual  loadings 
to  reservoirs.  In  accord  with  these  recommendations,  five  techniques 
for  estimating  the  annual  export  rate  of  total  phosphorus,  total  Kjeldahl 
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nitrogen,  and  suspended  solids  were  selected.  One  technique  was  based  on 
an  empirically  derived  equation  developed  for  the  U.  S.  Environmental 
Protection  Agency  (EPA)  by  Midwest  Research  Institute  (MRI)  (McElroy 
et  al.  1976).  A  regression  technique  developed  by  EPA's  National  Eutro¬ 
phication  Survey  (NES)  (Omernik  1977)  was  also  selected.  Additionally, 
three  methods  based  on  streamflow  and  referred  to  as  instream  techniques 
were  chosen:  average  flow,  average  concentration  (Q*C) ;  flow-weighted 
concentration  (QC);  and  flow  interval  (FI).  Evaluation  criteria  applied 
to  these  techniques  included  data  requirements,  necessity  of  local  exper 
tise,  and  relative  consistency  of  predictions.  Study  objectives  were 
to  determine  appropriate  applications  and  to  define  specific  advantages 
and  limitations  of  each  technique. 

7.  Since  one  of  the  primary  uses  of  these  techniques  may  be  to 
estimate  the  annual  loading  of  nutrients  to  reservoirs  for  subsequent 
use  with  simplified  techniques  for  predicting  trophic  status,  recommen¬ 
dations  were  included  for  the  application  of  each  technique  for  this 
purpose . 


PART  II:  ANNUAL  WATER  QUALITY  LOADING  TECHNIQUES 


8.  Detailed  descriptions  of  equations  comprising  each  of  the 
five  techniques  mentioned  in  paragraph  6  are  given  below.  Specific  data 
requirements,  coefficients,  and  references  are  also  included. 

Midwest  Research  Institute 

9.  The  objective  of  the  MRI  study  was  to  develop  procedures  for 
estimating  nonpoint  pollutant  sources  primarily  by  using  available  in¬ 
formation  and  local  technical  personnel. 

10.  The  MRI  Handbook  (McElroy  et  al.  1976)  includes  regionalized 
loading  functions  and  step-by-step  procedures  for  their  use.  In  addi¬ 
tion,  the  Handbook:  presents  graphically  and  in  tabular  form  some  of 
the  data  required  for  application  to  watersheds  throughout  the  United 
States;  provides  references  to  other  data  sources;  and  suggests  simple 
methods  for  generating  data  when  available  data  are  considered  inade¬ 
quate.  A  loading  function,  as  defined  by  the  MRI  Handbook,  is  a  "mathe¬ 
matical  expression  to  calculate  the  emission  of  a  pollutant  from  a  non¬ 
point  source  and  discharge  of  the  pollutant  into  surface  waterways." 

11.  The  basis  for  MRI  loading  functions  is  the  Universal  Soil 
Loss  Equation  (USLE),  developed  originally  by  W.  H.  Wischmeier  (Wisch- 
meier  and  Smith  1965)  to  describe  the  soil  transport  during  sheet  and 
rill  erosion  on  small  agricultural  fields.  The  USLE  expressed  the  amount 
of  soil  loss  based  on  the  soil  type,  slope  length,  crop  management,  and 
regional  rainfall  characteristics:* 

n 

Y(S)  =  T  A.  (RXKXLXSXCXPX  S  .) .  (1) 

i  d  i 

i=l 


*  See  the  MRI  Handbook,  which  contains  appropriate  figures  and  tables 
for  selecting  coefficients  and  necessary  conversion  factors  from  U.  S. 
customary  to  metric  (SI)  units. 
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where 


Y(S)  =  annual  yield  of  sediment  from  erosion  to  the  stream,  metric 
tons/yr 

n  =  number  of  subareas  in  the  watershed 
A^  =  acreage  of  watershed  subarea  i  ,  ha 

i  =  integer  from  1  to  n,  where  n  =  total  number  of  subareas 
i  in  a  watershed 

R  =  rainfall  erosivity  factor  (Figure  3-2,  MRI  Handbook),  ex¬ 
pressing  the  erosion  potential  of  average  annual  rainfall 
in  the  locality.  The  R  value  is  a  summation  of  the  in¬ 
dividual  storm  products  of  the  kinetic  energy  of  rainfall 
and  the  maximum  30-min  rainfall  intensity,  in  cm/hr,  for 
all  significant  storms,  on  an  average  annual  basis 

K  =  soil-erodibility  factor,  metric  tons/ha  for  each  R  unit 

L  =  slope-length  factor,  dimensionless 

S  =  slope-steepness  factor,  dimensionless 

C  =  cover  factor,  dimensionless 

P  =  erosion  control  practice  factor,  dimensionless 
=  sediment  delivery  ratio,  dimensionless 

12.  In  addition  to  sediment  load  from  erosion,  MRI  has  developed 
functions  and  data  requirements  that  depend  on  USLE  output  to  predict 
the  export  rate  of  total  nitrogen  and  total  phosphorus,  available  nitro¬ 
gen  and  phosphorus,  pesticides,  organic  matter,  heavy  metals,  traffic- 
related  pollutants,  and  acid-mine  drainage. 

13.  Total  nitrogen  load  has  two  components:  nitrogen  associated 
with  the  sediment  load  and  nitrogen  input  from  rainfall.*  Nitrogen  load 
can  be  determined  by: 


Y(N) 


Pr 


A  x 


Q(OR) 

Q(Pr) 


x  N. 


Pr 


x  b 


where 

Y(N)pr  =  stream  nitrogen  load  from  precipitation,  kg/yr 
A  =  area,  ha 


(2) 


*  See  the  MRI  Handbook,  which  contains  appropriate  figures  and  tables 
for  selecting  coefficients  and  necessary  conversion  factors  from  U.  S. 
customary  to  metric  (SI)  units. 
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Q(OR)  =  overland  flow  from  precipitation,  cm/yr 
Q(Pr)  =  total  amount  of  rainfall,  cm/yr 

Np^  =  nitrogen  load  in  precipitation,  kg/ha/yr 
b  =  attenuation  factor,  dimensionless 

The  total  nitrogen  loading  from  erosion  YCNT)^  is  estimated  by: 

Y(NT)  =  a  x  Y(S)  x  C  (NT)  x  r  (3) 

e  e  s  n 


where 

a  =  dimensional  constant 
YtS)^  =  total  sediment  load,  kg/yr 
C  (NT)  =  amount  of  nitrogen  in  the  top  30  cm  of  soil,  g/100  g 
&  =  0.55e-°08T  (1  -  e-°-005H) 

T  =  annual  average  temperature,  °C 
H  =  humidity  factor,  dimensionless 

and 


where 

Pr  =  precipitation,  cm/yr 

RH  =  relative  humidity,  percent 

SVP^  =  saturated  vapor  pressure  at  given  temperature, 
mm  of  Hg 

r  =  nitrogen  enrichment  ratio,  dimensionless 
n 

Total  phosphorus  load  is  estimated  by 

Y(PT)  =  a  x  Y(S)  x  C  (PT)  x  r  (5) 

e  s  p 


where 

Y(PT)  =  total  phosphorus  load  to  the  stream,  kg/yr 
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C  (PT)  =  amount  of  phosphorus  in  the  top  30  cm  of  soil,  g/100  g 
=  phosphorus  enrichment  factor,  dimensionless 

14.  The  MRI  Handbook  includes  appropriate  figures  and  tables  to 
determine  all  of  the  coefficients,  given  the  watershed  characteristics. 
These  characteristics  should  be  described  by  personnel  familiar  with  the 
watershed  soils,  land  use,  and  hydrology.  This  is  recommended  in  pref¬ 
erence  to  the  watershed  descriptions  provided  within  the  MRI  Handbook  or 
the  U.  S.  Geological  Survey  (USGS)  National  Atlas  (USGS  1970).  Geograph¬ 
ical  variability  between  watersheds  within  the  same  region  of  the  country 
is  not  represented  in  the  MRI  Handbook  and  National  Atlas  on  a  scale  com¬ 
mensurate  with  CE  planning  requirements  for  reservoirs.  The  MRI  Handbook 
is  based  on  data  from  first-order  watersheds,  whereas  the  National  Atlas 
describes  regional  areas,  e.g.,  major  river  basins. 

National  Eutrophication  Survey 

15.  The  EPA  conducted  a  nationwide  sampling  program  to  survey  this 
country's  lakes,  reservoirs,  and  tributaries  for  the  purpose  of  provid¬ 
ing  a  measure  of  their  trophic  status.  Based  on  the  sampling  of  more 
than  4000  NES  tributaries,  928  were  found  to  meet  the  nonpoint  source 
criterion  for  developing  a  nonpoint  source  loading  relationship.  Each 
watershed  was  sampled  once  a  month  for  a  year.  However,  the  sampling 
program  extended  over  a  3-year  period:  133  sites  were  sampled  in  the 
northeast  during  1972;  340  sites  were  sampled  in  the  east  and  southeast 
during  1973;  and  455  sites  were  sampled  in  the  west  and  midwest  during 
1974.  At  the  EPA  laboratory  in  Corvallis,  Oregon,  J.  M.  Omernik  (1977) 
developed  regional  regression  equations  for  lakes  and  reservoirs  by  re¬ 
lating  watershed  location,  land  use,  and  strearaflow  to  their  loading 
rates  of  phosphorus  and  nitrogen.  The  information  required  for  this 
method  included  the  general  location  of  the  watershed  and  the  percent 
land  use  of  the  watershed  classified  as  forest,  cropland,  and  urban. 

16.  Initially,  each  regression  equation  was  developed  using  char¬ 
acteristics  that  were  thought  to  influence  instream  nutrient  levels, 
e.g.,  soil  geology  and  land  use.  However,  it  was  found  that  land  use 
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was  the  watershed  characteristic  that  correlated  closer  to  stream  nutri¬ 
ent  concentrations.  Data  were  analyzed  with  respect  to  the  eastern, 
central,  and  western  regions  of  the  country.  Equations  are  listed  in 
Table  1  for  the  regional  regression  analyses  of  nitrogen  and  phosphorus 
and  related  interpretive  factors. 

17.  The  NES  regression  equation  estimates  a  mean  annual  concen¬ 
tration.  For  any  time  period  for  which  a  loading  value  is  required,  the 
stream  discharge  must  be  known  to  convert  the  mean  annual  concentration 
to  an  export  rate. 

18.  From  Table  1,  the  range  of  nutrient  concentrations  can  be 
determined  by  the  equation: 


Range  = 


_ f  x  r _ 

log^Q  nutrient  concentration 


(6) 


where 

f  =  multiplicative  standard  error 
r  =  correlation  coefficient 

log1Q  nutrient  concentration  =  value  estimated  by  equation 
The  range  (Equation  6)  was  found  to  contain  approximately  two  thirds  of 
the  observed  nutrient  concentrations  for  a  given  land  use  combination. 

A  first  approximation  of  the  mean  annual  concentration,  using  only  the 
land  use  information  for  a  selected  watershed,  can  be  obtained  from  the 
NES  regional  graphs  (Figures  23  and  24,  NES  Handbook).  In  addition, 
the  procedure  allows  for  refining  the  first  approximation  by  geograph¬ 
ical  area  when  nutrient  concentrations  are  expected  to  vary  from  the 
predicted  value  (Figures  25-28,  NES  Handbook  (Omernik  1977)).  A  more 
refined  estimate  of  the  mean  annual  nutrient  concentration  can  be  ob¬ 
tained  by  the  equation: 

Mean  annual  concentration  = 

(7) 

0  x  logjgf  +  logjQ  nutrient  concentration  (from  Table  1) 
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Table  1 

NES  Stream  Nutrient  Concentration  Predictive  Models  ( from  Omernik  1977) 


Region 

Model,  Correlation  Coefficient,  and  Multiplicative 
Standard  Error 

East 

L°810 

(PCONC)  = 

Total  Phosphorus 

1.8364  +  0.00971  (%  agric  + 
r  =  0.74,  f  =  1.85 

X  urb) 

Central 

Log10 

(PCONC)  = 

-1.5697  +  0.00811  (%  agric  • 
(X  for)  r  =  0.70,  f  =  2.05 

►  %  urb)  -0.002312 

West 

Log10 

(PCONC)  = 

1.1504  +  0.00460  (%  agric  + 
(%  for)  r  =  0.70,  f  =  1.91 

%  urb)  -0.00632 

Orthophosphorus 

East 

Log10 

(OPCONC)  = 

:  -2.2219  +  0.00934  (%  agric 
r  =  0.73,  f  =  1.86 

+  %  urb) 

Central 

Log10 

(OPCONC)  = 

:  -2.0815  +  0.00868  (%  agric 
r  =  0.63,  f  =  2.05 

+  %  urb) 

West 

Logio 

(OPCONC)  = 

:  -1.5513  +  0.00510  (%  agric 
(%  for)  r  =  0.64,  f  =  1.91 

+  %  urb)  -0.00476 

Total  Nitrogen 

East 

Log10 

(NCONC)  = 

-0.08557  +  0.00716  (%  agric 
(%  for)  r  =  0.85,  f  =  1.51 

+  %  urb)  -0.00227 

Central 

Logio 

(NCONC)  = 

-0.01609  +  0.00399  (%  agric 
(%  for)  r  =  0.77,  f  =  1.50 

+  %  urb)  -0.00306 

West 

Logio 

(NCONC)  = 

-0.03665  +  0.00425  (%  agric 
(%  for)  r  =  0.61,  f  =  1.75 

+  %  urb)  -0.00376 

Inorganic  Nitrogen 

East 

L°g10 

(INCONC)  = 

-0.3479  +  0.00858  (%  agric 
(X  for)  r  =  0.84,  f  =  1.93 

+  %  urb)  -0.00584 

Central 

Logio 

(INCONC)  = 

-0.5219  +  0.00482  (%  agric 
(%  for)  r  =  0.71,  f  =  2.06 

+  %  urb)  -0.00572 

West 

Log10 

(INCONC)  = 

-0.6339  +  0.00789  (%  agric 
(%  for)  r  =  0.65,  f  =  2.45 

+  X  urb)  -0.00657 
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where  a  =  standard  deviation  for  a  specific  geographic  area  (Figures  25 
27,  NES  Handbook) 

19.  An  annual  export  rate  for  a  specific  nutrient  can  be  esti¬ 
mated  as 


metric  tons/yr  =  C  *  31.536  *  F 


(8) 


where 

C  =  mean  annual  concentration,  mg/£ 

31.536  =  dimensional  constant 

F  =  mean  annual  flow,  cu  m/sec/yr 

20.  Omernik  (1977)  stated  that  the  relationships  indicate  corre¬ 
lations  only  between  existing  land  use  patterns  and  stream  nutrient  con¬ 
centrations.  These  factors  were  thought  to  explain  why  there  was  a  dif¬ 
ference  in  the  correlation  coefficient  for  each  of  the  models  describing 
nutrient  concent ra t ions  in  streams  from  the  three  regions  of  the  country 
Land  use  patterns  in  the  western  region  were  considered  significantly 
less  disturbed  by  man  than  the  eastern  region.  Hence,  the  correlation 
coefficient  was  higher  for  the  eastern  region  than  for  cither  of  the 
other  two  regions.  Predictions  of  loadings  using  these  regression  equa¬ 
tions  are  not  considered  valid  if  land  use  patterns  have  changed  signif¬ 
icantly  since  completion  of  the  NES. 

Average  Flow,  Average  Concentrat ion 

21.  This  approach  represented  perhaps  the  simplest  instream 
method  for  computing  annual  loadings  of  water  quality  constituents.  The 
mean  concentration  C  of  a  given  constituent  was  determined  using  all 
available  water  quality  data  i  collected  on  a  given  stream  or  river. 


n 


n 


(9) 
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where 


C  =  mean  concentration,  mg/fi 
i  =  an  integer  from  1  to  n 
n  =  number  of  measurements 

C  =  measured  concentration  for  each  data  measurement  i  ,  mg/£ 

1 

The  mean  flow  Q  was  determined  using  all  available  flow  data  i  on  a 
given  stream  or  river: 


n 


(10) 


where 

Q  =  mean  flow,  cu  m/sec 

Q.  =  flow  of  each  data  measurement  i  ,  cu  m/sec 
The  product  of  Equations  9  and  10,  Q*C  ,  represents  an  average  loading 
rate  over  a  year,  expressed  as  a  function  of  watershed  area.  This  is 
considered  to  be  an  estimate  of  the  annual  export  rate,  i.e.,  kilograms 
per  hectare  per  year.  The  data  required  for  this  method  are  usually 
available  in  the  CE  Districts  and  EPA's  STORET  (U.  S.  EPA  1979)  nation¬ 
wide  computer  library  of  water  quality  and  quantity  data  for  inland 
waters . 


Flow-Weighted  Concentration 


22.  Flow-weighted  loadings  QC  were  considered  essentially  as 

average  loading  rates.  Each  pair  of  flow  Q.  and  concentration  C. 

3  3 

data  is  multiplied  to  express  a  loading  rate  as  mass  per  unit  time. 
Next,  these  loading  rates  are  summed  over  the  data  set  and  divided  by 
the  number  of  values  n  : 


QC 


(11) 
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where 


j  =  each  data  pair 

n  =  number  of  values  within  the  data  set  for  which  flow  and 
concentration  values  have  been  recorded  simultaneously 

•=  product  of  flow  and  concentration,  mg/ cu  m  *  unit  time 

The  annual  export  rate  is  expressed  as  QC  summed  over  a  year  and  ex¬ 
pressed  as  a  function  of  the  watershed  area,  i.e.,  kilograms  per  hectare 
per  year. 

Flow  Interval 

23.  Initially,  the  flow  interval  (FI)  method  was  developed  in 
the  Lake  Erie  Wastewater  Management  Study  (U.  S.  Army  Engineer  District, 
Buffalo  1975)  to  estimate  the  phosphorus  loading  contributed  to  the  lake 
from  watersheds  with  little  or  no  water  quality  data.  Generally,  as  de¬ 
veloped,  the  method  should  be  considered  applicable  to  any  flow-dependent 
water  quality  constituent  of  interest.  It  requires  constituent  concen¬ 
tration  measurements  and  instantaneous  flow  data  (or  daily  flow  records) 
that  are  usually  available  from  the  USGS  District  Offices  or  CE  District 
Offices.  The  annual  water  quality  loading  is  determined  by  plotting  the 
product  of  a  constituent  concentration  times  the  flow  rate  at  the  time 

of  sampling,  i.e.,  discrete  flux  measurement,  versus  flow  for  the  period 

of  record.  Next,  the  number  of  flow  intervals  is  arbitrarily  determined 

at  n  and  the  flow  rate  Q.  at  the  end  of  each  interval  is  given  as 

1 

follows : 

Q 

Q  =  J5il* 
i  n 

where 

Q.  =  flow,  cu  m/sec,  at  the  end  of  each  interval 
i  =  flow  interval 

Q  =  largest  measured  flow,  cu  m/sec 
max 

n  =  number  of  flow  intervals 

24.  For  each  flow  interval  i  ,  the  mean  flux 
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m 


as  well  as  the  standard  error  of  the  mean 


2 
s  . 
1 


and 


k 

I 

2  i  =  l 


a.  -  l.) 


s  .  = 
l 


k(k  -  1) 


(13) 


(1A) 


where 

L.  =  discrete  flax  measurements  within  the  flow  interval  Q.  , 
mg/cu  m/sec 

k  =  number  of  flux  measurements 
The  mean  flux  L  for  the  period  of  record  is 


n 

L  =  S  L  P  (15) 

i  =  l 


and  the  variance  V  is 


V  = 


n 


I 

—  1 


s2P2 
1  1 


(16) 


where 


P.  =  probability  of  a  given  flow  occurring  in  flow  interval  i 
is  determined  from  the  available  flow  record  by 


P. 


(17) 


d  =  number  of  sampling  times  in  which  flows  in  interval  i 
occurred 

D  =  total  number  of  sampling  times  in  the  period  under 
consideration 
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The  range  of  values  for  the  average  flux  per  unit  time  L  for  a  given 
constituent  is 

L  =  L  ±  k  -jv~  (18) 

s 

where 

k  =  1.645  at  a  90  percent  confidence  interval  assuming  a  normal 
distribution 

For  an  annual  water  quality  export  rate  ,  the  average  flux  L  is 

multiplied  by  the  appropriate  time  period  t  to  convert  L  to  average 
flux  per  year: 

LA  =  t  x  L  (19) 

The  annual  water  quality  export  rate  represents  a  probabilistic 

flow-weighted  estimate  for  the  time  period  over  which  the  data  were 
collected. 
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PART  III:  SELECTION  OF  STUDY  WATERSHEDS 


25.  The  following  description  of  each  watershed  contains  the  in¬ 
formation  required  to  apply  and  evaluate  the  MRI  and  the  NES  regional 
loading  methods.  Instream  methods  only  require  water  quality  and  flow 
records;  however,  interpretation  of  the  loading  estimates  may  require 
watershed  information  similar  to  that  for  the  regional  methods. 

26.  Existing  data  sets  for  nonpoint  source  water  quality  load¬ 
ings  from  both  forested  and  agricultural  watersheds  were  used  for  eval¬ 
uating  the  selected  loading  techniques.  The  primary  criteria  were 
availability  of:  comprehensive  data  including  physiographic  features 
of  the  watershed;  flow  records;  and,  water  quality  parameters  of  sus¬ 
pended  solids,  total  phosphorus,  and  total  nitrogen. 

27.  The  selected  watersheds  represented  a  range  of  nonpoint 
sources  for  water  quality  constituents.  Included  were  the  forested 
Caddo  River  watershed  in  west  central  Arkansas  and  the  agricultural 
watersheds  of  the  Sandusky  River  and  Honey  Creek  within  the  Lake  Erie 
Basin  in  north  central  Ohio. 

Sandusky  River  Watershed* 

28.  The  Sandusky  River  is  approximately  185  km  long  and  drains 
approximately  368,000  ha  of  intensively  cultivated  agricultural  land  in 
north  central  Ohio  (Figure  1).  The  river  channel  shows  numerous  meanders 
and  flows  through  glacial  till  material.  The  channel  slope  ranges  from 
less  than  a  0.4-m/km  drop  to  a  4.7-m/km  drop  below  Tiffin,  Ohio.  Side- 
walls  of  the  channel  are  subject  to  sloughing  during  elevated  flow  con¬ 
ditions.  The  historical  flow  and  water  quality  records  (U.  S.  Army 
Engineer  District,  Buffalo  1978)  represent  an  extensive  data  collection 
effort  at  the  Fremont,  Ohio,  stream  gage. 

29.  The  Till  Plain  encompasses  the  southern  65  percent  of  the 


*  Descriptions  of  the  Sandusky  River  and  Honey  Creek  watersheds  were 
provided  by  the  U.  S.  Army  Engineer  District,  Buffalo  (1975). 
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Figure  1.  Sandusky  River  watershed  in  north  central  Ohio 
(from  U.  S.  Army  Engineer  District,  Buffalo  1975) 
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watershed  and  a  low  east-west  trending  moraine  crosses  the  watershed  just 
south  of  Tiffin,  Ohio.  The  northern  35  percent  of  the  watershed  is  lo¬ 
cated  within  the  Lake  Plains  subprovince.  The  surface  bedrock  formations 
consist  primarily  of  limestone,  dolomite,  and  shale.  The  presence  of 
groundwater  storage  is  small  since  the  formations  are  characteristically 
dense.  Glacial  drift  and  lake  deposits  are  thin  and  of  very  low  permea¬ 
bility.  Near  Bucyrus,  Crestline,  and  Upper  Sandusky,  the  glacial  drift 
is  highly  permeable  and  thicker  in  the  vicinity  of  the  moraines. 

30.  Glacial  deposits  resulted  in  the  development  of  very  complex 
soil  types  with  a  variable  texture.  In  the  northern  35  percent  of  the 
watershed,  the  soils  are  predominantly  the  very  poorly  drained  Hoytville 
silty  clay  loam  and  clay.  In  the  southern  65  percent  of  the  watershed, 
the  soils  are  moderately  fine-textured  Wisconsin  glacial  till.  The 
western  part  of  the  Till  Plains  consists  of  soils  developed  over  lime¬ 
stone,  whereas  the  soils  in  the  eastern  portion  are  underlain  by  sand¬ 
stone  and  shale.  The  groundwater  contribution  to  streamflow  is  quite 
small  with  a  low  sustained  flow.  However,  the  agricultural  land  in  the 
western  portion  of  the  watershed  has  been  made  productive  only  through 
provision  of  underground  tile  drainage  to  allow  more  productive  farming. 

31.  Very  detailed  land  use,  soils,  and  topographic  information 
were  provided  through  computer  analysis  of  high  altitude  aerial  photo¬ 
graphic  imagery.  This  information  was  developed  as  part  of  the  Buffalo 
District's  Lake  Erie  Wastewater  Management  Study  (U.  S.  Army  Engineer 
District,  Buffalo  1975).  A  subsequent  project  (Cahill,  Pierson,  and 
Cohen  1978),  concerning  evaluation  of  best  management  practices,  as 
demonstrated  on  the  Honey  Creek  watershed,  contributed  to  this  very  de¬ 
tailed  data  set  describing  land  use  and  soils.  The  study  area  repre¬ 
sented  about  3,200  sq  km  above  the  city  of  Fremont,  Ohio,  and  the  USGS 
gage.  Of  the  approximately  320,000  ha  included  in  this  study,  256,200  ha 
(80  percent)  is  comprised  of  cropland;  7,300  ha  (2  percent)  is  pasture; 
28,700  ha  (9  percent)  is  urban;  21,400  ha  (7  percent)  is  forest;  and  the 
remaining  6,400  ha  (2  percent)  is  surface  water.  Major  crops  include 
corn,  soybeans,  and  wheat.  Approximately  200,000  people  reside  in  this 
watershed . 
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Honey  Creek  Watershed* 


32.  Honey  Creek  is  a  tributary  to  the  Sandusky  River  (Figure  2). 
Approximately  48,400  ha  is  drained  by  this  creek.  The  channel  is  72  km 
long  and  slopes  an  average  of  3  m/km.  The  existing  data  file  reflects 
39,100  ha  of  watershed  above  the  Melmore,  Ohio,  stream  gage.  The  his¬ 
torical  flow  and  water  quality  data  collection  was  initiated  in  1976. 

The  data  include  numerous  elevated  flow  sampling  and  water  quality  anal¬ 
yses  of  samples  obtained  every  4  hr. 

33.  Soil  mapping  of  this  watershed  was  completed  by  the  U.  S. 
Department  of  Agriculture  (USDA)  Soil  Conservation  Service  (SCS)  and  the 
Ohio  Department  of  Natural  Resources.  The  soils  are  fine  textured  al¬ 
though  quite  variable,  e.g.,  muck  to  sand.  The  watershed  lies  in  the 
glacial  till  region  of  the  Sandusky  River  watershed;  hence,  the  soils 
have  low  permeability  and  very  poor  natural  drainage.  The  seasonal 
water  table  is  within  15  cm  of  the  land  surface  for  approximately  82  per¬ 
cent  of  the  area.  Approximately  80  percent  of  the  land  has  less  than 

2  percent  slope  and  only  4  percent  of  the  land  has  a  slope  greater  than 
6  percent.  An  important  characteristic  of  the  soils  is  their  inherently 
high  erodibility  coupled  with  long  slope  lengths  to  a  channel. 

34.  Tile  drainage  systems  are  in  use  throughout  the  watershed  to 
increase  field  drainage  and  enable  mechanized  fieldwork  earlier  in  the 
spring  and  later  in  the  fall. 

35.  Of  the  39,100  ha  total  area,  31,800  ha  (81  percent)  is  crop¬ 
land.  Major  crops  include  corn,  soybeans,  and  wheat.  Approximately 
4,200  ha  (11  percent)  is  forested,  2,000  ha  (5  percent)  is  urban,  800  ha 
(2  percent)  is  wetlands,  and  300  ha  (1  percent)  is  water. 

Caddo  River  Watershed 


36.  The  Caddo  River  watershed  is  in  west  central  Arkansas 


*  The  description  of  the  Honey  Creek  watershed  was  provided  by  the 
U.  S.  Army  Engineer  District,  Buffalo  (1975). 
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RELATION  TO  COUNTIES 
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Figure  2.  Honey  Creek  watershed  in  north  central  Ohio 
(from  U.  S.  Army  Engineer  District,  Buffalo  1975) 
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(Figure  3).  The  river  flows  a  distance  of  126  km  from  its  source  in  the 
Ouachita  Mountains  to  its  confluence  with  the  Ouachita  River.  The 
watershed  has  an  area  of  117,600  ha  with  an  elevation  range  from  675  m 
in  the  mountains  on  the  west  to  53  m  in  the  east  at  the  confluence  with 
the  Ouachita  River. 


•  MT  IDA 


Figure  3.  Caddo  River,  Arkansas,  watershed  showing 
Glenwood  gage 

37.  The  headwaters  region  of  about  52,000  ha  lies  to  the  west- 
northwest  of  Glenwood,  Arkansas,  and  is  part  of  the  Novaculite  Uplift. 
This  region  is  composed  of  long,  eastward-trending,  even-crested  moun¬ 
tains  and  flat  intermountain  basins.  The  region  is  underlain  by  mix¬ 
tures  of  shale  and  sandstone  rocks.  However,  some  chert  is  observable  on 
the  south  side  of  the  watershed  with  a  few  outcrops  on  the  north  side. 

38.  West  of  Glenwood,  Arkansas,  the  dominant  soil  associations 
are  the  Carnasaw,  Pickens,  and  Pickwick.  The  Pickens  soils  are  found 
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on  steep  mountains  and  ridges.  They  are  characterized  as  very  shallow, 
acid,  excessively  drained,  and  moderately  permeable  loamy  soils  with 
considerable  gravel.  The  Carnasaw  soils  are  found  on  rolling  hills  and 
are  moderately  deep,  acid,  and  well  drained.  They  are  derived  from  a 
quartzite  and  sandstone  bedrock  with  an  overlying  clay  loam  subsoil. 

The  Pickwick  soils  are  found  near  the  streambanks  and  are  typically 
deep,  acid,  and  well  drained  with  a  lot  of  gravel  near  the  surface. 

39.  Forested  land  occupies  more  than  80,000  ha  (68  percent)  of 
the  watershed  with  oak,  hickory,  and  pine  being  the  dominant  species. 
Agriculture  occupies  about  35,300  ha  (30  percent)  of  the  watershed  and 
consists  of  pastures,  hay  meadows,  and  abandoned  fields  in  varying 
stages  of  succession.  The  remaining  2  percent  of  the  watershed  area 
is  urban  and  developed  land  with  a  population  estimated  to  be  less 
than  1500  people. 

40.  The  climate  of  the  watershed  is  generally  mild.  Summers  are 
quite  long  with  temperatures  occasionally  greater  than  40°C,  and  winters 
are  usually  short  and  moderate.  Average  annual  precipitation  is  135  cm 
with  43  percent  of  the  precipitation  resulting  in  runoff. 

41.  Runoff  tends  to  be  higher  in  the  headwater  portion  of  the 
watershed  west  of  Glenwood,  Arkansas.  This  portion  of  the  watershed, 
which  comprises  only  44  percent  of  the  total  area,  contributes  over 
60  percent  of  streamflow.  The  stream  gage  at  Glenwood,  Arkansas,  was 
selected  for  this  study  since  it  had  the  best  historical  flow  record 
(26  years)  and  several  years  of  water  quality  data. 

42.  The  peak  discharge  frequency  curve  for  the  Caddo  River  at  the 
Glenwood  gage  is  shown  in  Figure  4.  The  mean  annual  peak  storm  flow, 
interpolated  from  this  curve,  is  about  500  cu  m/sec.  The  range  of  an¬ 
nual  peak  flows  during  the  period  of  record,  1940-1976,  was  from 

2500  cu  m/sec  in  1968  to  80  cu  m/sec  in  1941,  which  corresponds  to 
water  depths  of  approximately  10  m  and  3  m,  respectively,  at  the  Glen¬ 
wood  gage. 


25 


PART  IV:  APPLICATION  AND  COMPARISON  OF  TECHNIQUES  FOR  COMPUTING 
NONPOINT  SOURCE  ANNUAL  WATER  QUALITY  EXPORT  RATES 


Introduction 

43.  The  data  requirements,  ease  of  application,  and  assumptions 
and  limitations  characteristic  of  each  technique  are  described  in  this 
section.  Techniques  were  categorized  as  follows:  (a)  those  regionalized 
techniques  requiring  physiographic  information  about  the  watershed,  and 
(b)  those  techniques  that  require  instream  water  quality  and  flow  mea¬ 
surements.  Regionalized  techniques  included  MRI  (McElroy  et  al.  1976) 
and  NES  (Omernik  1977).  Instream  techniques  included  the  average  flow, 
average  concentration  (Q*C);  flow-weighted  concentration  (QC);  and  flow 
interval  (FI).  Comparison  of  these  techniques  using  identical  data  sets 
provides  the  information  required  to  objectively  select  the  technique(s) 
most  suitable  for  a  specific  application. 

44.  Comparison  and  evaluation  of  results  within  the  regionalized 
techniques  category  emphasized  the  following:  (a)  quantity  and  detail 
of  land  use,  soil,  and  meteorological  information  required;  (b)  optimum 
watershed  size  for  its  application;  (c)  type  of  watershed  for  best  ap¬ 
plication,  i.e.,  forested  or  agricultural;  and  (d)  assumptions  and  limi¬ 
tations  of  the  techniques.  These  topics  were  addressed  by  comparing 
effects  of  the  physiographic  information  requirements  on  the  estimated 
annual  water  quality  export  rate  for  the  Caddo  River  watershed.  The 
first  evaluation  used  only  general  physiographic  information  found  in 
the  Phase  I  General  Design  Memorandum  (U.  S.  Army  Engineer  District, 
Vicksburg  1966)  for  the  DeGray  Lake;  the  second  evaluation  included  very 
detailed  land  use,  soil,  and  topographic  analyses  that  had  been  completed 
from  a  previous  investigation.*  The  assumptions  and  limitations 

*  J.  Nix  et  al.  1974.  "Collection  of  Environmental  and  Water  Quality 
Data  on  DeGray  Lake  and  Caddo  River  Watershed"  (unpublished).  Annual 
Contract  Report,  prepared  by  Arkansas  Water  Resources  Research  Center, 
Fayetteville,  Ark.,  for  U.  S.  Army  Engineer  Waterways  Experiment  Sta¬ 
tion,  CE,  Vicksburg,  Miss.,  under  Contract  No.  DACW39-73-C-0125 ;  re¬ 
ports  prepared  in  1975,  1976,  and  1977  under  Contract  No.  DACW39-75- 
C-0025 . 
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associated  with  the  NES  and  MRI  techniques  were  discussed  to  determine 
the  applicability  of  each  technique  to  forested  and  agricultural 
watersheds . 

45.  The  instream  techniques  category  was  evaluated  based  on  ade¬ 
quacy  of  historical  flow  and  water  quality  records,  sampling  frequency, 
and  advantages  and  limitations  associated  with  their  application.  These 
techniques  were  used  to  estimate  annual  export  rates,  using  available 
flow  and  water  quality  data  for  each  watershed. 

Midwest  Research  Institute 

46.  The  MRI  Handbook  provided  a  step-by-step  procedure  with  guid¬ 
ance  on  specific  data  requirements  for  determining  mean  annual  export 
rates  of  a  variety  of  pollutants  from  nonpoint  sources.  The  MRI  proce¬ 
dures  evaluated  in  this  study  included:  annual  soil  loss  by  erosion  to 
waterways,  total  phosphorus  (TP)  annual  export  rate,  and  total  nitrogen 
annual  export  rate.  For  this  evaluation,  total  Kjeldahl  nitrogen  (TKN) 
was  defined  as  the  MRI's  total  nitrogen  minus  nitrate-nitrogen  and 
nitrite-nitrogen  values.  As  previously  stated,  all  export  rates  should 
be  considered  as  long-term  average  annual  loadings  to  the  stream  during 
a  time  period  without  significant  land  use  changes. 

Caddo  River  watershed 

47.  Generalized  physiographic  information  (Table  2),  as  described 
in  the  Phase  I  General  Design  Memorandum  (U.  S.  Army  Engineer  District, 
Vicksburg  1966)  for  DeGray  Lake  and  MRI  Handbook  information,  was  used 
in  the  first  evaluation  to  determine  mean  annual  export  rates  for  soil, 
TP,  and  TKN.  Table  3  was  developed  from  appropriate  tables  and  figures 
in  the  MRI  Handbook. 

48.  The  most  difficult  items  to  estimate  without  a  thorough  anal¬ 
ysis  of  the  topography  were  the  overall  percent  slope  of  the  land  area 
and  the  slope  length.  The  slope  length  was  defined  in  the  MRI  Handbook 
as  "the  average  distance  for  overland  flow  prior  to  intersecting  a  stream 
channel."  No  guidance  was  provided  in  the  MRI  Handbook  to  determine 
slope  length.  Therefore,  10  first-order  watersheds  in  agricultural  and 
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Table  2 

General  Description  of  Caddo  River  Watershed 


_ Item _ 

Watershed  area,  % 
Area ,  ha 

Slope  of  land,  % 
Slope  length,  m 
General  information 


_ Agriculture _ 

20 

10,500 

6 

61 

Grass  and  shrub  cover 
over  approx  80%  of 
area 

Light  to  moderate 
grazing 


_ Forest _ 

80 

42,000 

13 

46 

Tree  canopy,  70% 

Litter  cover,  80% 

Managed  undergrowth 

Medium  stocked  with 
timber 


MR  I 

Table  3 

Coefficient  Values  Specific  to 

Caddo  River  Watershed* 

Land  Use  Areas 

Coefficient 

Agriculture 

Forest 

A 

10,500 

42,000 

LS 

0.9 

2.6 

C 

0.013 

0.004 

Sd 

0.15 

0.35 

DD 

0.07 

0.37 

*  From  U.  S.  EPA  (1979). 
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forested  areas  were  arbitrarily  selected  on  a  topographic  map  of  the 
Caddo  River  watershed  to  represent  the  minimum  and  maximum  slope  length 
and  percent  slope.  The  length  and  percent  slope  were  determined  for  the 
10  watersheds  and  averaged  to  represent  the  overall  watershed  area. 

49.  The  land  use  cover  factor  C  and  sediment  delivery  ratio  S. 

d 

were  subjectively  estimated.  To  obtain  an  accurate  estimate  of  these 
factors  would  require  more  detailed  information  than  normally  would  be 
available.  An  estimate  of  C  may  be  obtained  from  the  appropriate 
tables  in  the  MRI  Handbook  with  only  a  general  knowledge  of  the  watershed 
land  use.  Coefficient  is  determined  after  estimating  the  drainage 

density  DD  ,  which  is  the  ratio  of  stream  length  to  watershed  area 
within  each  land  use  classification.  For  the  Caddo  River  watershed, 

2 

it  was  estimated  that  approximately  70  stream  kilometres  drained  104  km 

2 

of  agriculture  and  96  stream  kilometres  drained  420  km  of  forest.  The 
S^  was  difficult  to  estimate  because  of  its  ambiguous  association  with 
DD  (Figure  3-10,  MRI  Handbook).  It  is  assumed  that  the  rate  of  sediment 
delivery  from  the  land  surface  to  a  stream  is  dependent  primarily  on  the 
number  of  streams  per  unit  watershed  area. 

50.  A  more  detailed  description  of  the  Caddo  River  watershed  is 
given  in  Table  4.  This  information  was  obtained  from  comprehensive  land 
use  maps  and  soils  classification  data  provided  under  contract  by  the 
Arkansas  Water  Resources  Research  Center*  and  by  the  Arkansas  Soil  Con¬ 
servation  Service.  Inis  detailed  physiographic  information,  along  with 
the  appropriate  figures  and  tables  in  the  MRI  Handbook  and  the  National 
Atlas  (USGS  1970),  allowed  a  more  thorough  evaluation  and  subsequent  se¬ 
lection  of  coefficient  values  representative  of  specific  soil  types  and 
land  uses  (Tables  5  and  6).  Subsequent  determination  of  mean  annual  ex¬ 
port  rates  for  soil,  TKN,  and  TP  provided  a  comparison  between  the  use  of 
the  generalized  and  detailed  data  to  estimate  the  annual  export  rate  of 


Nix  et  al . ,  op .  cit . 
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Table  4 


Detailed  Description  of  Caddo  River  Watershed 
Above  Glenwood,  Arkansas,  Gage 


Range- 

Agri- 

Cleared 

Item 

Forest 

land 

culture 

Unproductive 

Urban 

Watershed  area,  % 

86 

6 

3 

4 

1 

Area ,  ha 

45,150 

3150 

1575 

2100 

525 

Slope  of  area,  % 

12 

3 

3 

3 

2 

Slope  length,  m 

40 

50 

60 

50 

60 

General  information: 

Canopy  cover,  % 

70 

30 

25 

35 

Ground  cover,  % 

85 

50 

30 

35 

Table  5 

Coefficient 

Values*  for  Determining  Annual  : 

Sediment  Yield 

from  Caddo 

River  Watershed 

Land  Use-Soil  Type 

A 

K 

LS 

C 

Sd 

Forest 

Carnasaw 

18,060 

0.46 

0.70 

0.003 

0.78 

Sherwood 

9,030 

0.26 

0.40 

0.003 

0.52 

Pickens 

18,060 

0.48 

1.50 

0.002 

0.56 

Rangeland 

Ca  rnasaw 

1,575 

0.46 

0.35 

0.008 

0.78 

Sherwood 

1,575 

0.26 

0.25 

0.008 

0.52 

Agriculture 

Carnasaw 

315 

0.46 

0.18 

0.200 

0.78 

Pickwick 

1,260 

0.46 

0.14 

0.100 

0.78 

Cleared  unproductive 

Sherwood 

1,260 

0.26 

0.35 

0.013 

0.52 

Pickens 

840 

0.48 

0.35 

0.13 

0.56 

Urban 

Carnasaw 

105 

0.46 

0.30 

0.010 

0.78 

Sherwood 

367 

0.26 

0.35 

0.010 

0.52 

Pickwick 

53 

0.46 

0.20 

0.010 

0.78 

*  From  U.  S.  EPA  (1979)  and  USGS  (1970). 


Table  6 


Constants*  Used  to  Determine  Water  Quality  Loadings 
from  Caddo  River  Watershed 


Constant 

Value 

R 

330 

P 

1.0 

Q(0R) 

0.15 

Q(Pr) 

1.5 

Pr 

b 

1.0 

For  estimating  water  quality 

Temp,  °C 

15.8 

loadings  from  precipitation 

H 

348 

Pr ,  mm 

1408 

RH,  % 

70 

SVPt 

13.48 

r 

2.0 

n 

r 

1.5 

p 

Cg(NT),  g/100  g 

0.13 

Cg(PT),  g/100  g 

0.02 

*  From  U.  S.  EPA  (1979)  and  USGS  (1970). 


soil,  TKN,  and  TP  from  the  Caddo  River  watershed: 


General  physiographic  data, 
kg/ha/yr  (metric  tons/yr) 

Detailed  physiographi c  data, 
kg/ha/yr  (metric  tons/yr) 

Reduction,  percent 


Soil 

723(39,929) 

622(32,675) 

14 


_ TKN  _ 

1.95(102) 

1.59(84) 

18 


TP 

0.23(12) 

0.19(9.8) 

17 


A  change  of  approximately  14  to  18  percent  in  the  estimated  annual  export 
rate  of  soil,  TKN,  and  TP  for  the  Caddo  River  watershed  was  obtained  us¬ 
ing  the  detailed  physiographic  data.  Hence,  the  development  of  detailed 
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physiographic  data  may  not  be  justified  in  some  cases,  based  on  this  one 
application. 

51.  Since  the  MRI  method  is  based  on  the  relation  of  water  quality 
loadings  to  soil  loss  and  subsequent  suspended  solids  (SS)  production, 
coefficient  values  for  determining  sediment  yield  must  be  selected  with 
caution.  In  addition,  the  nitrogen  and  phosphorus  concentrations  in  the 
top  30  cm  of  soil,  C  (NT)  and  C  (PT)  ,  must  be  considered  critical 
parameters.  An  error  of  one  to  two  orders  of  magnitude  can  be  intro¬ 
duced,  depending  on  the  selected  concentration.*  Man's  activities,  e.g. 
cultivated  fields,  managed  forests,  and  urban  areas,  may  contribute  sig¬ 
nificantly  to  nitrogen  and  phosphorus  export  from  a  specific  watershed; 
therefore,  Cs(NT)  and  Cs(PT)  should  be  carefully  evaluated.  It  is 
recommended  that  local  SCS  personnel  familiar  with  a  specific  watershed 
be  contacted  for  assistance  in  making  appropriate  estimates. 

Sandusky  River  and 
Honey  Creek  watersheds 

52.  Watershed  characteristics  required  for  application  of  the  MRI 
method  are  given  in  Tables  7  and  8.  Descriptions  of  watershed  physiog¬ 
raphy  were  provided  by  the  Buffalo  District.  This  information  was  de¬ 
veloped  from  analyses  of  high  altitude  photographs  and  topographic  maps 
supplemented  by  ground  truth  data  collection  and  evaluation. 

53.  Based  on  available  information,  the  Honey  Creek  watershed  was 

considered  representative  of  the  land  use,  relief,  and  soils  of  the  San¬ 
dusky  River  watershed.  Likewise,  the  percent  slope  and  slope  length  for 
the  Sandusky  River  and  Honey  Creek  watersheds  were  considered  the  same 
for  this  comparison.  Only  the  and  A  were  different  (Table  9). 

The  remaining  coefficient  values  that  were  selected  for  the  Sandusky 
River  and  Honey  Creek  watersheds  are  listed  in  Table  9.  Only  loading  by 
erosion  was  considered  in  the  computations.  Rainfall  data  for  these 
watersheds  were  not  available  from  the  Buffalo  District  nor  the  U.  S. 
Weather  Bureau. 


*  The  USGS  based  these  soil  nutrient  concentrations  on  regional  geo- 
morphological  soil  development  patterns. 
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Table  7 

Description  of  Sandusky  River  Watershed 


_ I tem _ 

Watershed  area,  %* 
Area ,  ha 

Slope  of  area,  % 
Slope  length,  m 


Cropland 

82 

263,408 

2 

172 


Woodland 

9 

28,749 

4 

307 


General  information 


Corn,  soybeans,  and  wheat  Tree  canopy,  70% 
cover  80%  of  the 

watershed  area  Litter  cover,  80% 


Some  erosion  control 
management 

Ground  cover,  50% 


Managed  undergrowth 

Medium  stocked  with 
timber 


*  The  remaining  9  percent  represents  other  land  use  activities  and 
water . 


Table  8 

Description  of  Honey  Creek  Watershed 


Item 

Cropland 

Woodland 

Watershed  area,  %* 

82 

10 

Area ,  ha 

32,052 

4096 

Slope  of  area,  % 

2 

4 

Slope  length,  m 

172 

307 

General  information 

Corn,  soybean,  and  wheat 

Tree  canopy,  70% 

cover  51%;  alfalfa 
clover  and  fescue 

Litter  cover,  80% 

cover  31%  of  the 

Managed  undergrowth 

watershed  area 

Medium  stocked  with 

Some  erosion  control 
management 

Ground  cover,  60% 


*  The  remaining  8  percent  represents  other  land  use  activities  and 
water . 


Table  9 


MR1  Coefficients  Used  to  Determine  Annual  SS ,  TKN , 
and  TP  Export  Rates  from  Sandusky  River 
and  Honey  Creek  Watersheds 


MR  I 

Sandusky  River 

Watershed 

Honey  Creek  Watershed 

Coefficient 

Cropland 

Woodland 

Cropland 

Woodland 

A 

263,408 

28,749 

30,052 

4096 

K 

0.33 

0.33 

0.33 

0.33 

LS 

0.44 

0.77 

0.40 

0.70 

C 

0.60 

0.003 

0.65 

0.003 

Sd 

0.21 

0.17 

0.28 

0.22 

R 

120 

120 

120 

120 

P 

0.37 

1.0 

0.37 

1.0 

r 

n 

3.0 

2.0 

3.5 

2.0 

r 

p 

1.5 

1.5 

1.5 

1.5 

C  (NT) 

0.30 

0.30 

0.18 

0.18 

C  (PT) 
s 

0.15 

0.15 

0.15 

0.15 

54.  The  mean 

i  annual  export 

rates  of  SS 

,  TKN,  and  TP  for  the  San- 

dusky  River  and  Honey  Creek  watersheds  were  calculated  using 

the  selected 

coefficient  values 

in  Table  9.  Results  are  as  follows: 

Mean  Annual  Export 

Rates,  kg/ha/y 

r  (metric 

tons/yr ) 

Watershed 

SS 

TKN 

TP 

Sandusky  River 

1824  (480 

1 , 553) 

11.0  (2883) 

4.1  (1081) 

Honey  Creek 

1958  (76, 

807) 

10.6  (415) 

4.4  (173) 

The  annual  loading 

estimates  varied  primarily 

as  a  function 

of  watershed 

area.  Apparently, 

small  changes  in  a  few  coefficients,  i.e. 

Sd  '  rn  ’ 

r  ,  C  (NT)  ,  and 
P  s 

C  (PT)  ,  have 
s 

negligible 

impact  on  the  MR1  results. 

Discussion 

55.  Since  the  coefficients 

are  multiplicative  in  the 

MRI  equa- 

tions,  any  error  in 

i  specific  coefficients  is 

multiplied  to  a 

larger  error 
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Two  tables  in  the  MR 1  Handbook,  Tables  3-10  and  4-3,  provide  a  range  of 
loading  values  for  si spended  solids,  nutrients,  and  organic  matter.  The 
extreme  values  of  each  coefficient,  based  on  available  field  plot  data, 
were  used  by  MRI  to  calculate  the  probable  range  shown  in  Tables  3-10  and 
4-3  of  that  report.  These  tables  illustrate  the  best  estimate  of  the 
annual  "loading  range  of  accuracy”  using  this  method  and  available  sup¬ 
porting  data.  Based  on  the  loading  estimates  for  the  Caddo  River,  San¬ 
dusky  River,  and  Honey  Creek  watersheds,  the  probable  range  of  accuracy 
in  export  estimates  is  presented  as  follows: 

Range  in  Export  Estimates,  kg/ha/yr  (metric  tons/yr ) 
Watershed  _ Soil _  TP  _ TKN 

Caddo  River  10-1,000  0. 1-3.0  0.1-10 

(52-52,550)  (5-158)  (5-525) 

Sandusky  River  100-5,000  2-10  5-20 

(26,341-1,317,055)  (527-2634)  (1317-5628) 

Honey  Creek  100-5,000  2-10  5-20 

(39,200-195,550)  (78-392)  (196-784) 

The  range  of  loadings  for  soil,  TP,  and  TKN  is  so  broad  as  to  warrant 
omitting  the  MRI  procedure  by  itself  for  estimating  loadings  to  reser¬ 
voirs  by  CE  Districts. 

56.  Advantages .  A  summary  of  the  advantages  of  the  MRI  Handbook 
method  that  are  pertinent  to  CE  planning  functions  is  presented  as 
follows : 

a.  Requires  minimal  field  data. 

b.  Provides  opportunity  to  solicit  local  expert  assistance 
for  selecting  appropriate  coefficient  values. 

c.  Provides  for  the  use  of  all  available  information  on  a 
specific  watershed  for  considering  appropriate  coefficient 
values . 

d.  Requires  minimal  manpower,  time,  and  cost  to  provide  esti¬ 
mates  of  loadings  for  input  into  early  project  planning 
and  design. 

e.  May  provide  insight  into  effects  of  changing  land  use  by 
changes  in  coefficient  values,  e.g.,  C  . 

f.  Is  based  on  the  USLE,  which  is  a  primary  analytical  tool 
extensively  used  throughout  the  United  States  by 
agriculturists . 
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57.  Limitations .  A  summary  of  the  limitations  of  the  MRI  method 


follows : 

a.  The  USLE  was  developed  for  field-size  plots,  not  large 
watersheds.  The  suitable  aerial  scale  for  application  has 
not  been  determined. 

b.  The  selection  of  USLE  coefficient  values  is  highly 
subjective. 

c.  All  chemical  loadings  are  associated  with  sediment  yield 
(suspended  solids).  This  is  highly  dependent  on  the  nu¬ 
trient,  mineral,  and  organic  matter  content  of  watershed 
soils . 

d.  The  USLE  is  comprised  of  terms  lacking  theoretical  basis, 

e.g.,  sediment  delivery  ratio  and  slope  length,  which 
cannot  be  measured  in  the  field  for  large  watersheds. 

e.  The  temporal  variation  of  C  is  not  documented.  Seasonal 
variations  in  watershed  cover  are  not  discussed  or  incor¬ 
porated  in  the  MRI  technique. 

f.  No  explicit  consideration  is  made  in  the  USLE  for  runoff 
quantity.  Terms  implicitly  relating  to  rainfall  are  R  , 

K  ,  LS  ,  and  C 

g.  The  seasonal  effect  of  ra inf a  1 1/ runof f  quantity  is  not 
considered.  The  R  values  represent  the  portion  c f  the 
year  when  rainfall  intensity  and  frequency  is  highest  and 
evapotranspiration  is  lowest.  For  example,  in  northern 
Ohio,  80  percent  of  the  erosive  rainfall  occurs  from  May 
through  October,  yet  most  of  the  runoff  occurs  from  Novem¬ 
ber  through  April  (U.  S.  Army  Engineer  District,  Buffalo 
1975).  In  Arkansas,  runoff-producing  storms  occur  primar¬ 
ily  from  January  through  May  (U.  S.  EPA  1979).  As  defined, 
most  water  quality  export  occurs  during  these  months  and 
not  throughout  the  year,  as  MRI ' s  estimated  annual  export 
rate  suggests. 

h.  A  variation  of  at  least  10  to  100  tiroes  can  be  introduced 
in  the  annual  loading  estimates  of  total  nitrogen  and 
total  phosphorus  if  C  (NT)  and  C  (PT)  are  in  error 
for  watershed  subarea  i  . 


National  Eutrophication  Survey 


58.  The  NES  technique  consists  of  regionalized  regression  equa¬ 
tions  incorporating  three  general  categories  of  land  use:  forest,  agri¬ 
culture,  and  urban.  The  rationale  and  data  used  to  develop  this  nonpoint 
source  loading  technique  are  thoroughly  described  in  the  NES  Handbook 
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(Omernik  1977).  Unlike  the  MRI  methodology,  consideration  of  annual 
sediment  production,  i.e.,  suspended  solids,  is  not  included  in  this 
procedure . 

59.  The  NES  regional  equations  were  used  to  estimate  the  mean 
annual  concentrations  of  orthophosphorus  and  total  phosphorus  and  inor 
ganic  and  total  nitrogen  (TN)  in  the  Caddo  River,  Sandusky  River,  and 
Honey  Creek.  The  only  data  required  were  the  percentage  of  each  land 
use  (Table  10),  and  mean  annual  flow  for  each  stream. 

Table  10 

Percentage  Land  Use  Composition  of  Caddo  River,  Sandusicy  River, 
and  Honey  Creek  Watersheds 


Watershed 

Agriculture 

Forest 

Urban 

Caddo  River 

3 

86 

1 

Sandusky  River 

82 

9 

6 

Honey  Creek 

82 

10 

5 

60.  The  estimated  standard  error  a  for  water  quality  data  rep¬ 
resentative  of  the  Caddo  River,  Sandusky  River,  and  Honey  Creek  water¬ 
sheds  was  selected  from  Figures  24-27  in  the  NES  Handbook.  The  suggested 
range  of  standard  error  (Table  11)  for  each  geographical  area  was  also 
used  to  compare  the  respective  range  of  estimated  loading  values  for  the 
water  quality  parameters.  The  selection  of  a  standard  error  is  subjec¬ 
tive;  it  is  based  on  the  assumption  that  a  given  watershed  would  be  sim¬ 
ilar  to  the  average  watershed  for  that  specific  geographical  region.  In 
developing  the  regional  regression  equations,  the  NES  team  assumed  that 
they  were  sampling  from  a  normal  distribution  of  watersheds  without 
point  source  pollution.  However,  only  928  of  the  4000  NES  tributary 
sites  met  this  criterion.  The  broad  range  of  standard  error,  as  speci¬ 
fied  for  each  geographical  region  of  the  United  States,  is  illustrative 
of  the  water  quality  variability  found  within  each  region. 

61.  The  mean  annual  nutrient  concentrations  (Table  12)  were  sim¬ 
ilar  to  the  NES  results,  presented  in  Figures  3  and  4  of  the  NES  Hand¬ 
book,  for  the  estimated  TP  concentration  in  >75  percent  forested  and 
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Table  11 

Estimated  Standard  Error  a  (mg/&)  and  Suggested  Range*  of 
Selected  Water  Quality  Parameters  for  Caddo  River,  Sandusky 
River,  and  Honey  Creek 


Water  Quality 

Caddo  River 

Sandusky  River 

Honey  Creek 

Parameter 

a 

Range 

a 

Range 

a 

Range 

Total 

phosphorus 

0.1 

-0.5  to  0.5 

0.1 

-0.5 

to  0.5 

0.1 

-0.5 

to  0.5 

Ortho¬ 

phosphorus 

0.3 

-0.5  to  0.5 

-0.4 

-0.5 

to  0.5 

-0.4 

-0.5 

to  0.5 

Total  nitrogen 

-0.2 

-0.5  to  0.5 

1.0 

0.5 

to  1.0 

0.8 

0.5 

to  1.0 

Inorganic 

nitrogen 

-0.5 

-1.0  to  -0.5 

0.2 

-0.5 

to  0.5 

0.1 

-0.5 

to  0.5 

*  NES  Handbook  (Omernik  1977). 


Table  12 

Comparison  of  Annual  TP  and  TKN  Concentrations 
with  NES  Regression  Analyses 


Watershed 


Caddo  River 

Sandusky  River 

Honey  Creek 

NES  (>75% 
forest) 

NES  (>75% 

agriculture) 


TP,  mg /SL 


TKN ,*  mg/£ 


Estimate 

0.02 

0.12 

0.11 


Range 


(0.01-0.03) 

(0.07-0.12) 

(0.07-0.14) 


Estimate 

0.41 

2.3 

2.5 


Range 


(0.36-0.59) 
(1.8-4. 7) 
(1.0-4. 3) 


0.02 


0.14 


0.52 


1.13 


*  TKN  =  TN  -  (N03-N  +  N02-N). 
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agricultural  watersheds.  Good  agreement  was  also  observed  for  TKN 
(equivalent)  in  the  forested  watersheds;  however,  the  estimated  TKN 
(equivalent)  in  the  agricultural  watersheds  was  twice  that  estimated  by 
NES  analysis  for  similar  watersheds.  Two  reasons  for  the  overestima¬ 
tion  may  be  considered:  (a)  the  seasonal  loading  patterns  of  nitrogen 
cannot  be  examined  as  described  in  the  NES  Handbook  since  the  regional 
relationships  were  developed  on  only  a  1-year  data  set;  and  (b) 
nitrogen  loading  depends  significantly  on  time  of  year,  quantity  and 
method  of  fertilizer  applications,  and  tilling  operations. 

62.  Maps,  presented  in  Figures  25-28  of  the  NES  Handbook,  were 
developed  to  assist  in  selecting  an  appropriate  standard  error  represen 
tative  of  the  study  regions.  However,  they  were  not  sufficiently  de¬ 
tailed  to  represent  all  basins  within  a  region.  Only  in  those  regions 
where  NES  data  were  concentrated  was  it  possible  to  estimate  a  reason¬ 
able  standard  error.  This  may  explain  why  the  range  of  loadings  for 
phosphorus  and  nitrogen  as  estimated  for  the  Caddo  River  was  so  close 
to  the  regional  loading  means;  whereas,  the  range  of  loadings  for  the 
Honey  Creek  and  Sandusky  River  watersheds  was  not  close  to  the  regional 
loading  means. 

63.  Table  13  shows  the  estimated  average  annual  export  rate  of 
TP  and  TKN  for  the  Caddo  River,  Sandusky  River,  and  Honey  Creek  water¬ 
sheds.  Significantly  higher  nutrient  loadings  are  estimated  from  the 


Table  13 

Comparison  of  Annual  Water  Quality  Export  Rates  Using  NES  Regional  Analysis* ** 


kg/ha/yr 

Range 


kg/ha/yr 

Range 


etric  tons 


kg/ha/yr  metric  tons/yr  kg/ha/yr  metric  tons 

Watershed  Avg  Range  Avg  Range  Avg  Range  Avg  Rap8 

Caddo  River  0.17  (0.11-0.23)  9  (5-14)  3.5  (3. 1-5.1)  186  (163-268) 

Sandusky  River  0.77  (0.45-0.77)  188  (110-188)  15  (12-30)  3605  (2821-7366) 

Honey  Creek  0.48  (0.33-0.62)  19  (12-24)  11  (4-19)  434  (173-746) 


*  Range  assuming  recommended  lower  and  upper  limits  of  standard  error.  The  NES 
computed  ranges  are  listed  in  parentheses,  i.e.,  based  on  ±1  standard  error. 

**  TKN  =  TN  -  (N03-N  +  N02-N). 
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agricultural  watersheds  than  from  the  forested  watersheds.  The  computed 
range  (in  parentheses)  was  based  on  the  NES's  standard  error  for  the 
regions . 

64.  In  comparison  to  the  MRI  analysis  (paragraph  56),  a  slightly 
lower  average  and  a  narrower  range  in  annual  TP  export  estimates  were 
computed  using  the  NES  technique.  However,  the  annual  TKN  export  esti¬ 
mates  for  the  MRI  and  NES  methods  were  quite  similar.  From  these  re¬ 
sults,  both  techniques  compared  favorably  for  the  forested  Caddo  River 
watershed  only. 

65.  It  would  appear  that  the  NES  analysis  provided  a  better  plan¬ 
ning  tool  than  the  MRI  method  for  estimating  loadings  because  it  is 
based  on  actual  data  from  watersheds  throughout  the  United  States.  How¬ 
ever,  the  NES  study  specified  that  the  regional  regression  equations 
were  applicable  only  to  the  data  used  in  their  analyses  and  the  spe¬ 
cific  time  covered  by  the  respective  sampling  period.  Since  the  data 
were  not  flow-weighted  and  represent  only  1  year  of  monthly  data  collec¬ 
tion  for  each  region,  the  mean  annual  nutrient  concentration  expressed 
for  each  watershed  in  the  study  does  not  relate  to  year-to-year,  month- 
to-month,  or  storm-to-storm  variability  (McDowell  and  Omernik  1979). 
Subsequent  changes  in  land  use  since  completion  of  the  NES  sampling 
program  could  greatly  alter  the  coefficient  values  and  land  use  factor 
composition  within  each  regionalized  equation.  Finally,  there  was  no 
empirical  basis  for  the  NES  loading  relationships  as  was  purported  by 
the  MRI  method,  which  uses  the  USLE.  The  NES  procedure  would  appear  to 
be  most  reliable  for  large,  geographically  similar  watersheds  as  shown 

in  Figures  25-28  of  the  NES  Handbook.  For  these  reasons  and  those  listed 
below,  CE  application  of  this  loading  technique  for  planning  is  not  rec¬ 
ommended  except  when  applied  to  watersheds  in  geographical  areas  where 
a  large  number  of  NES  watersheds  exist. 

66.  Advantages .  Advantages  of  the  NES  method  are  the  following: 

a.  Only  percentage  land  use  composition  of  a  watershed  and 
mean  annual  flow  are  required  to  apply  the  technique. 

b.  Required  figures  and  tables  are  provided  in  the  NES 
Handbook. 
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c.  Provision  is  made  for  computing  a  range  in  mean  annual 
loadings  of  total  nitrogen  and  inorganic  nitrogen. 

d.  Very  little  cost  and  effort  are  required  to  obtain 
loading  values. 

e.  There  is  no  limit  to  the  size  of  watershed  for  applica¬ 
tion  of  this  loading  technique. 

67.  Limitations .  The  limitations  of  the  NES  method  are  as 
follows : 

a.  All  variability  between  watersheds  within  a  region  is  ex¬ 
pressed  as  percentage  land  use  differences.  Furthermore, 
heavy  reliance  is  placed  on  available  aerial  photography 
and  maps  to  delineate  land  uses.  Many  sites  do  not  have 
this  information  available. 

b.  Agricultural  and  urban  land  uses  are  not  differentiated 
relative  to  nonpoint  source  loadings  of  nutrients. 

c.  A  historical  flow  record  is  required  to  obtain  a  realis¬ 
tic  mean  annual  flow  without  being  biased  toward  wet-  or 
dry-year  periods. 

d.  Regionalized  equations  may  have  shown  correlations  be¬ 
tween  existing  land  use  and  water  quality  concentration 
patterns  during  the  NES  study.  No  theoretical  basis 
exists  for  predicting  loadings  if  a  change  in  land  use  or 
water  quality  concentration  occurred  since  the  NES  com¬ 
pleted  their  sampling  program. 

e.  Large  areas  of  the  United  States  were  not  included  in  the 
study,  e.g.,  southwest  and  coastal  plains  of  the  southeast. 

f .  Only  nitrogen  and  phosphorus  loadings  were  included  in 
the  NES  analysis. 


Q*C  ,  QC  ,  and  FI  Instream  Techniques 

Q*C  technique 

68.  Instream  techniques  have  been  used  extensively  within  CE  Dis¬ 
tricts  to  compute  annual  nonpoint  source  water  quality  loadings  to  res¬ 
ervoirs.  The  Q*C  method  has  been  found  to  be  predominant  based  upon 
a  telephone  survey  of  each  CE  District.  Since  most  instream  data  col¬ 
lection  programs  have  been  conducted  in  cooperation  with  USGS  and  have 
utilized  USGS  procedures,  the  majority  of  the  available  data  reflect  the 
prevalent  flow  and  stream  water  quality  conditions  associated  with  base 
flow.  Moreover,  surface  water  quality  data  for  these  streams  have  been 
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based  on  weekly,  biweekly,  or  monthly  grab  samples.  Water  quality  data 
representative  of  elevated  flow  conditions,  i.e.,  incorporating  hydro¬ 
graph  sampling,  are  not  adequately  considered  in  most  routine  sampling 
programs  to  determine  the  influence  of  elevated  flow  on  constituent  ex¬ 
port  rates.  Consequently,  suspended  solids  and  nutrient  loadings  based 
on  historical  data  are  biased  toward  low-flow  conditions.  Similar  con¬ 
clusions  were  made  by  a  previous  investigation  (Johnson  et  al.  1976). 

This  technique  is  most  appropriate  for  estimating  annual  water  quality 
loadings  of  conservative  constituents ,  e.g.,  pH,  sodium,  and  potassium, 
where  only  data  from  routine  sampling  programs  are  available. 

QC  technique 

69.  The  QC  technique  has  been  used  for  determining  the  annual 
export  rate  of  both  flow-dependent  and  flow- independent  water  quality 
constituents.  However,  application  of  this  technique  does  not  improve 
the  export  rate  estimates  for  those  constituents ,  e.g.,  TKN,  associated 
with  organic  components  of  runoff  unless  the  data  are  representative  of 
elevated  and  base  flow  conditions.  The  composition  and  quantity  of  par¬ 
ticulate  and  dissolved  organic  components  in  runoff  are  influenced  by 
many  factors,  including  season;  storm  intensity,  duration,  and  frequency; 
and  degree  of  decomposition  of  the  organic  material  (Lystrom  et  al.  1978, 
Johnson  et  al.  1976,  Fisher  et  al.  1968,  Biggar  and  Corey  1969).  The 
relative  influence  of  these  factors  on  the  computed  annual  loading  esti¬ 
mates  depends  on  the  number  of  samples  and  period  of  time  over  which  the 
samples  are  collected  and  analyzed.  To  lessen  the  influence  of  these 
factors,  sampling  is  required  to  adequately  represent  the  seasonal  dif¬ 
ferences  and  characteristic  storm  patterns  in  a  specific  watershed. 
Consequently,  the  design  of  the  sampling  program  must  consider  the  fre¬ 
quency  of  occurrence  of  elevated  flows  and  associated  water  quality 
loadings  prior  to  applying  this  technique. 

FI  technique 

70.  The  FI  technique  (Fisher  et  al.  1968)  is  essentially  a  stream 
water  quality  rating  method.  It  is  based  on  identifying  the  frequency 
of  occurrence  of  different  flows  and  determining  the  water  quality  load¬ 
ing  rate  for  specified  flow  intervals.  The  technique  uses  instantaneous 
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flow  and  water  quality  data  for  computing  the  flow-weighted,  mean  load¬ 
ing  rates  within  selected  flow  intervals.  The  impact  of  infrequent 
high-  and  low-flow  conditions  and  associated  constituent  concentrations 
on  the  estimated  water  quality  loading  rate  is  reduced  significantly  by 
weighting  the  computed  export  rate  by  the  frequency  of  occurrence  of  each 
flow  value  in  the  historical  flow  record.  Though  the  nature  and  amount 
of  field  data  required  for  application  will  vary  from  project  to  project, 
elevated  and  low  flow  conditions  must  be  represented.  Consequently,  his¬ 
torical  data  and/or  additional  field  sampling  must  include  samples  ob¬ 
tained  over  a  wide  range  of  streamflow.  A  cost-effective  instream  sam¬ 
pling  program  should  be  designed  to  sample  throughout  selected  flow 
hydrographs  during  different  seasons  of  the  year.  Guidance  on  sampling 
frequencies  is  provided  by  Johnson  (1979).  Data  should  include  either 
discrete  concentration  measurements  and  daily  flow  records  or  discrete 
concentration  measurements  and  simultaneous  flow  measurements.  Appro¬ 
priate  application  of  the  instream  techniques  depends  on:  (a)  the 
quality  and  extent  of  the  historical  flow  and  water  quality  data, 

(b)  the  adequate  representation  in  the  data  set  of  low  and  elevated 
streamflow  conditions,  and  (c)  the  amount  of  funds  and  time  available 
to  obtain  field  data,  if  sufficient  data  do  not  exist  for  specific  proj¬ 
ects.  Where  seasonal  changes  may  influence  the  export  rate  of  selected 
chemical  constituents,  e.g.,  orthophosphate  and  nitrite-  and  nitrate- 
nitrogen,  a  better  estimate  of  the  annual  export  rate  is  obtained  by 
grouping  the  data  and  computing  seasonal  export  rates. 

Discussion 

71.  The  influence  of  elevated  versus  low  flows  on  the  estimated 
loading  rates  for  the  Caddo  River  is  shown  in  Table  14.  Approximately 
200  samples  from  nine  storm  flow  hydrographs  on  the  Caddo  River  were 
analyzed.  The  storms  represented  winter,  summer,  and  fall  seasons. 

High-  and  low-flow  data  were  separated  arbitrarily  into  two  groups  for 
analysis,  i.e.,  data  above  or  below  the  26-year  mean  daily  flow  of  !'♦  cu 
m/sec.  Approximately  70  percent  of  the  200  samples  were  representative 
of  flows  less  than  14  cu  m/sec.  The  flow  frequency,  as  characterized  by 
nine  flow  bydrographs,  corresponded  to  the  26-year  flow  record  for  the 


**  No.  of  samples:  67  elevated  flow,  107  low  flow,  174  all  flow 
t  No.  of  samples:  67  elevated  flow,  157  low  flow,  224  all  flow 
tt  One  standard  error  expressed  as  percent. 


Glenwood  gage.  Approximately  80  percent  of  the  measured  flows,  or 
7680  samples,  were  below  14  cu  m/sec. 

72.  The  standard  errors  of  the  estimates  were  computed,  assuming 
that  statistically  different  loading  relationships  existed  for  SS ,  TKN, 
and  TP  at  flows  less  than  or  greater  than  14  cu  m/sec.  Only  the  Q*C 
and  QC  loading  techniques  were  compared,  since  the  FI  technique,  by 
definition,  is  not  applicable  to  this  type  of  comparison.  The  impor¬ 
tance  of  elevated  flows  in  exporting  water  quality  constituents  was 
clearly  illustrated  by  comparing  the  export  rates  and  standard  error 
(Table  14) .  The  highest  loading  of  SS  ,  TKN  ,  and  TP  occurred  duri.ig 
flows  >14  cu  m/sec.  Moreover,  the  SS  export  rate  computed  using  the  QC 
technique  was  approximately  three  times  the  export  rate  computed  using 
the  Q*C  method.  Only  after  combining  the  low  and  elevated  flow  data 
were  the  estimated  loading  and  percent  standard  errors  significantly  im¬ 
proved  for  both  methods.  These  results  indicate  a  need  for  a  represen¬ 
tative  data  set  incorporating  low  and  elevated  flow  conditions  at  their 
statistically  normal  frequency  of  occurrence,  if  the  Q*C  and  QC  tech 
niques  are  to  be  used. 

73.  The  estimated  mean  annual  export  rates  and  standard  error  of 
the  estimates  for  SS,  TKN,  and  TP  for  the  Caddo  River,  Sandusky,  and 
Honey  Creek  watershed?-  (Figure  5)  were  compared.  The  number  of  flow  and 
water  quality  samples  used  for  this  analysis  is  listed  in  Table  15.  The 
mean  annual  water  quality  export  rates  computed  by  the  Q*C  and  QC 
techniques  varied  within  only  a  factor  of  two  to  three  for  each  water¬ 
shed  (Figure  5).  However,  when  the  water  quality  export  rates  were 

Table  15 

Number  of  Samples  Analyzed  Using  the 
Instream  Loading  Techniques 


Watersheds 

Flow 

Suspended 

Solids 

Total  Kjeldahl 
Nitrogen 

Total 

Phosphorus 

Caddo  River 

226 

195 

174 

224 

Sandusky  River 

1337 

1253 

194 

1289 

Honey  Creek 

1202 

1181 

276 

1191 

Figure  S.  Comparison  of  estimated  mean  annual  export 
rates  using  mstream  data 

flow-weighted  and  adjusted  for  their  probability  of  occurrence  (FI 
technique),  the  estimated  annual  export  rate  and  standard  error  of  est i 
mate  were  reduced  significantly.  Moreover,  the  influence  1  infrequent 
elevated  flows  and  subsequent  high  water  quality  constituent  concentra¬ 
tions  on  computed  annual  export  rates  (Figure  5,  FI  and  FI)  was  r  -dm 

)  2. 

significantly  compared  to  the  QC  technique. 

74.  The  FI  technique  yielded  a  twolold  reduction  in  the  computed 
annual  export  rate  for  SS ,  TKN ,  and  TP  as  compared  with  the  other  in- 
stream  loading  methods  (Table  14  and  Figure  5).  Also,  the  standard 
error  of  the  estimate  was  reduced  to  19  to  21  percent  for  SS  and  8  to 
11  percent  for  TKN  and  TP.  Similar  results  were  reported  for  the  tri¬ 
butaries  of  Lake  Erie  by  the  U.  S.  Army  Engineer  District,  Buffalo 
(1975).  Furthermore,  there  was  virtually  no  change  in  the  computed 
annual  export  rate  and  error  of  the  estimate  by  using  a  2-cu-m/sec  tl.w 
interval  rather  than  the  5-cu-m/sec  interval  (Table  14  and  Figure  v 

75.  The  export  rates  of  SS ,  TKN,  and  TP  from  the  Sandusky  Rivet 
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and  Honey  Creek  watersheds  (Figure  5)  were  nearly  identical.  Estimated 
SS  and  TP  export  rates  using  the  FI  method  were  similar  to  the  respec¬ 
tive  export  rates  computed  by  the  Q*C  and  QC  methods.  An  adequate 
number  of  samples  (Table  15)  had  been  analyzed  for  these  rivers  to 
statistically  represent  the  natural  frequency  of  occurrence  for  low- 
and  e leva ted- f 1 ow  conditions.  The  TKN  export  rate,  as  determined  by  the 
FI  method,  was  much  higher  for  Honey  Creek  than  for  the  Sandusky  River. 
Moreover,  the  export  rate  of  SS  from  the  Sandusky  River  watershed  was 
more  than  double  the  rate  from  the  Honey  Creek  watershed.  Explanations 
for  these  differences  may  include  any  number  of  factors:  unique  point 
sources,  land  use  and  land  use  management  differences  in  the  watersheds, 
soil  type  differences,  and  areal  variation  of  rainfall  in  the  watersheds. 
In  constrast,  the  annual  export  rates  of  TP  from  the  two  agricultural 
watersheds  were  nearly  identical. 

76.  The  export  rates  of  TP,  SS ,  and  TKN  using  the  FI  analysis 
(Figure  5)  were  similar  to  previous  research  comparing  runoff  from  for¬ 
ested  and  agricultural  watersheds  (True  1976,  Fisher  et  al.  2968,  and 
Timmons,  Burwell,  and  Hoit  1973).  The  estimated  annual  export  rate  of 
TP  was  10  times  lower  for  the  Caddo  River  watershed  than  for  either  the 
Sandusky  River  or  Honey  Creek  watershed.  The  SS  export  rate  from  the 
Honey  Creek  watershed  was  twice  the  rate  from  the  Caddo  River  watershed. 
The  TKN  export  from  the  forested  Caddo  River  watershed  was  twice  the 
rate  of  either  the  agricultural  Sandusky  River  or  Honey  Creek  watersheds. 

77.  The  appropriate  application  of  the  instream  techniques  for 
estimating  mean  annual  export  rates  depends  on:  the  quality  and  extent 
of  the  historical  flow  and  water  quality  data;  the  adequate  representa¬ 
tion  in  the  data  set  of  low  id  elevated  streamflow  conditions;  and  the 
amount  of  time  and  funds  available  to  obtain  additional  data  for  spe¬ 
cif  i  <.  projects.  The  data  used  for  evaluating  these  loading  techniques 
represented  very  extensive  sampling  programs.  One  objective  of  the  sam¬ 
pling  program  for  each  watershed  was  to  obtain  sufficient  i ns t ant aneous 
flow  and  water  quality  data  to  permit  water  quality  constituent  loading 
analyses.  Unfortunately,  most  of  the  flow  records  from  gaging  stations 
throughout  t  fie  United  States  report  mean  daily  flows,  whereas  water 


quality  data  represent  instantaneous  measurements.  The  estimated  mean 
annual  export  rate  using  available  data  should  not  bias  the  estimate  to 
an  unacceptable  level  as  long  as  high  flows  are  represented  (U.  S.  Army 
Engineer  District,  Buffalo  1975).  However,  this  estimate  could  be  mis¬ 
leading  only  if  the  mean  daily  flows  significantly  underestimate  water 
quality  loadings  occurring  at  intervals  more  frequent  than  daily,  i.e., 
during  storm  events  (Westerdahl,  Perrier,  and  Nix  1977). 

78.  It  is  difficult  to  specify  the  adequacy  of  a  data  set.  The 
respective  elevated-  and  low-flow  conditions  represented  by  the  data 
should  be  consistent  with  the  intended  application.  In  the  FI  analysis, 
it  was  assumed  that  storm  events  were  most  important  in  the  loading  of 
water  quality  constituents  from  the  watershed.  For  example,  the  approxi¬ 
mately  200  samples  (Table  15)  analyzed  for  the  Caddo  River  watershed 
included  9  storm  hydrographs  that  represent  winter,  summer,  and  fall 
storms.  Since  most  of  the  SS ,  TKN ,  and  TP  loadings  occurred  at  elevated 
flow  (Table  19),  the  data  were  considered  representative  of  the  water 
quality  loadings  from  all  watershed  sources. 

79.  These  instream  techniques  are  not  appropriate  for  computing 
loading  estimates  for  periods  other  than  the  period  of  record  being 
analyzed.  However,  the  FI  method  was  considered  the  most  reasonable 
technique  based  on  advantages  and  limitations  listed  in  Table  16.  Gen¬ 
erally,  any  instream  technique  that  incorporates  the  probability  of 
occurrence  for  specific  flows  and  water  quality  concentrations  would  be 
most  appropriate  for  estimating  mean  annual  water  quality  export  rates 
for  the  time  period  associated  with  the  data. 

General  Discussion 

80.  Mean  annual  loadings  of  SS,  TKN,  and  TP  for  the  Caddo  River, 
Sandusky  River,  and  Honey  Creek  watersheds  are  summarized  in  Table  17. 

The  MRI  method  provides  the  broadest  range  of  loadings  for  SS,  TKN,  and 
TP.  Selection  of  each  MRI  coefficient  involves  some  degree  of  subjec¬ 
tivity.  Any  errors  made  in  each  coefficient  selection  are  multiplied 
during  computations  in  solving  the  soil  loss  equation.  Moreover,  TP 
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Table  16 


Characteristic  Advantages  and  Limitations 
of  Inst  ream  Techniques 


_  Advantages _ _  Linn  tat  i  on s_ 

Average  Flow,  Average  Concentration 

Mean  annual  flow  and  water  quality  con-  Water  quality  constituent  export  rate 
stituent  concentrations  are  computed  is  assumed  to  be  flow  independent, 

using  all  available  historical  data 
for  the  stream,  e.g.,  mean  daily  flows 
and  discrete  water  quality  measurements. 

Minimal  manpower,  time,  and  cost  are  re¬ 
quired  to  provide  estimates  of  loading 
for  input  into  early  project  planning 
and  design. 


Flow-Weighted  Concentration 


Only  the  historical  flow  data  are  re¬ 
quired  at  the  time  samples  for  water 
quality  were  analyzed.  Any  addi¬ 
tional  field  studies  would  be  limited 
to  sampling  at  different  flow  rates 
to  obtain  these  data. 

Minimal  manpower,  time,  aod  cost  are 
required  to  provide  estimates  of 
loadings  for  input  into  early 
project  planning  and  design. 


Water  quality  constituents  are  assumed 
to  be  flow  dependent. 

Frequency  of  elevated  and  low  flows  is 
not  a  consideration.  Consequently, 
the  mean  annual  loading  rate  may 
either  be  higher  or  lower  than  the 
actual  loading  rate  for  the  stream. 


Flow  Interval 


This  is  valid  for  most  flow-dependent 
and  f low- independent  water  quality 
constituents . 

Limited  data  that  represent  storm 

hydrograph  sampling  characteristic  of 
different  storm  intensities  and  dura¬ 
tion  throughout  a  year  are  adequate 
for  computing  realistic  annual  water 
quality  loadings  from  a  watershed. 


Data  bases  are  required  that  include 
either  discrete  concentration  mea¬ 
surements  and  the  daily  flow  record 
or  discrete  concentration  measure¬ 
ments  and  simultaneous  flow 
measurements . 

Considerable  data  representative  ol 
seasonal  loadings  are  required  tor 
those  constituents  influenced  by 
seasonal  changes,  e.g.,  orthophos¬ 
phate  and  nitrite-  and  nitrate- 
nitrogen.  The  annual  export  rate  of 
these  constituents  can  be  determined 
by  grouping  the  data  by  season  and 
using  the  FI  method  to  compute  an 
export  rate  for  each  season. 

Though  selection  of  flow  interval  size 
is  arbitrary,  considerable  effort  may 
be  required  to  group  the  data  accord¬ 
ingly.  This  is  especially  significant 
if  the  data  base  is  not  computerized. 
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Table  17 


Comparison  of  Mean  Annual  Loadings  of  SS ,  TKN,  and  TP  for  the 
Caddo  River,  Sandusky  River,  and  Honey  Creek  Watersheds 


Mean  Annual  Loading,  metric  tons/yr 


Loading  Technique 

SS 

TKN 

TP 

Caddo  River 

MRI  (No . 1 ) 

39,929 

102 

12 

(No. 2) 

32,675 

84 

10 

Range 

52-52,550 

5-158 

5-525 

NES 

-- 

215  (+24%) 

10  (+50%) 

Q*C 

18,770  (+23%) 

418  (+18%) 

19  (+19%) 

QC 

35,448  (+30%) 

536  (+17%) 

34  (+24%) 

FI  @  5  cu  m/sec 

16,036  (+21%) 

280  (+9%) 

19  (+8%) 

FI  @  2  cu  m/sec 

13,523  (+20%) 

263  (+7%) 

17  (+8%) 

Sandusky  River 

MRI 

480,553 

2883 

1081 

Range 

26,341-1,317,055 

1317-5628 

527-2634 

NES 

-- 

7527  (+11%) 

150  (+30%) 

Q*C 

209,675  (+9%) 

2240  (+11%) 

464  (+7%) 

QC 

396,097  (+8%) 

3979  (+11%) 

749  (+7%) 

FI  @  5  cu  m/sec 

337,782  (+4%) 

520  (+3%) 

632  (+2%) 

FI  @  2  cu  m/sec 

273,420  (+5%) 

395  (+3%) 

506  (+2%) 

Honey  Creek 

MRI 

76,807 

415 

173 

Range 

39,200-195,550 

196-784 

78-392 

NES 

-- 

492  (+62%) 

18  (+33%) 

Q*C 

13,923  (+10%) 

353  (+9%) 

45  (+8%) 

QC 

34,721  (+17%) 

666  (+9%) 

83  (+8%) 

FI  @  5  cu  m/sec 

26,408  (+9%) 

142  (+8%) 

67  (+6%) 

FI  @  2  cu  m/sec 

19,907  (+12%) 

142  (+7%) 

57  (+7%) 

and  TKN  export  rates  are  subject  to  the  error  of  the  SS  loading  estimate 
as  well  as  the  range  of  values  for  C^CPT)  and  C^CNT).  Most  impor¬ 
tantly,  however,  the  MRI  technique  requires  that  the  user  be  familiar  with 
watershed  land  use,  soil,  and  hydrologic  characteristics.  The  MRI  pro¬ 
cedure  does  offer  a  technique  for  estimating  changes  in  water  quality 
loadings  as  a  function  of  future  changes  in  land  use. 

81.  The  NES  method  of  regional  regression  analysis  has  only 
limited  use  for  CE  application.  Only  those  CE  project  areas  within  the 
general  vicinity  of  the  NES  study  sites  can  be  considered  appropriate. 
Loading  estimates  should  be  obtained  by  selecting  the  most  appropriate 
standard  deviation  represented  in  Figures  24-28  of  the  NES  Handbook  for 
a  specific  watershed. 

82.  Of  the  three  instream  techniques,  the  FI  method  is  recom¬ 
mended.  Though  historical  flow  and  water  quality  records  are  desired, 
analysis  of  even  limited  seasonal  data  can  be  performed  with  the  tech¬ 
nique.  Interval  size  has  only  a  slight  effect  on  the  computed  annual 
loadings  and  standard  error  of  the  estimate  (Table  17). 
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PART  V:  IMPLICATIONS  FOR  ASSESSING  RESERVOIR  EUTROPHICATION 


Introduction 

83.  Throughout  project  planning  and  during  the  operation  of  exist¬ 
ing  projects,  recurring  questions  are  raised  concerning  the  trophic  sta¬ 
tus  of  a  reservoir.  Vollenweider  (1968)  proposed  one  of  the  first  meth¬ 
ods  based  on  a  correlation  of  total  annual  point  and  nonpoint  source 
loadings  versus  lake  response  rather  than  on  internal  lake  processes. 
However,  this  method  was  developed  using  a  data  base  primarily  from 
northern,  temperate  natural  lakes  that,  therefore,  is  not  applicable 

for  most  CE  reservoirs.  Several  investigators,  e.g.  Shannon  and  Brezonik 
(1972),  Dillon  (1974,  1975),  Dillon  and  Rigler  (1974),  and  Larsen  and 
Mercier  (1976),  attempted  to  develop  better  techniques  or  to  modify 
Vollenweider ' s  methodology.  These  methods,  also  based  on  the  use  of 
external  nutrient  loadings,  i.e.,  phosphorus,  provide  an  evaluation  of 
the  eutrophication  state  of  existing  lakes. 

84.  The  TP  loadings,  estimated  by  each  technique  in  this  study, 
were  used  to  illustrate  what  could  happen  when  the  loading  estimates  are 
used  to  assess  reservoir  eutrophication.  The  eutrophication  assessment 
method  proposed  by  Larsen  and  Mercier  (1976)  was  selected  to  compare  the 
effects  of  various  annual  loading  estimates.  This  method  uses  mean  an¬ 
nual  phosphorus  loading,  net  annual  phosphorus  sedimentation,  and  general 
reservoir  morphometric  and  hydraulic  considerations  to  provide  an  esti¬ 
mate  of  trophic  status,  e.g.,  oligotrophic ,  mesotrophic,  and  eutrophic. 
However,  many  assumptions  and  limitations  must  be  recognized  in  order 
that  this  or  similar  procedures  may  be  used  and  interpreted  appropri¬ 
ately.  Some  of  these  were  summarized  by  Larsen  and  Mercier  (1976)  and 
Dillon  (1974).  The  various  methods  make  the  following  assumptions: 

a.  Phosphorus  concentration  is  an  indication  of  lake  trophic 
state.  Phosphorus  is  considered  to  be  the  limiting 
nutrient . 

b.  Rate  of  water  inflow  and  outflow  is  steady  and  equal. 
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c.  The  lake  is  an  instantaneously  well-mixed  system;  i.e., 
inflows  are  instantaneously  dispersed  throughout  the 
lake . 

d.  There  is  a  net  loss  of  phosphorus  to  the  sediment  annually. 

e.  Sedimentation  rates  are  at  steady-state  in  the  lake. 

f.  Outflow  water  constituent  concentrations  are  representa¬ 
tive  of  the  mean  concentrations  in  the  lake. 

g.  With  a  change  in  loading  rate,  the  lake  biomass  will 
increase  or  decrease  accordingly. 


Application  and  Evaluation  of  Annual  Loading  Techniques 
for  Assessing  Reservoir  Trophic  State 


85.  Of  the  three  watersheds  studied,  only  the  data  set  for  the 
Caddo  River  contained  the  requisite  data.  Therefore,  it  was  selected 
for  the  comparative  evaluation. 

86.  To  estimate  the  trophic  status  of  DeGray  Lake,  the  Larsen 
and  Mercier  (1976)  relation  was  used: 


P  = 


1  -  R 


exp 


(20) 


phosphorus  concentration  in  the  reservoir,  mg/£ 

incoming  phosphorus  concentration  at  steady-state,  mg/£ 

empirically  determined  phosphorus  retention  coefficient, 
i.e., 

R  =  0.854  -  0.142  In  q  (21) 

exp  s 

q^  =  areal  hydraulic  load,  dimensionless,  i.e, 

q  =  Z  x  a  (22) 

s  w 

Z  =  mean  reservoir  depth,  m  (13.3  m  for  DeGray  Lake) 

o  =  theoretical  hydraulic  washout  coefficient, 
w 


where 


P  = 
P  = 


exp 
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The  theoretical  hydraulic  washout  coefficient  is: 


a 

w 


T 


where 


=  hydraulic  residence  time,  i.e., 

^  _  reservoir  volume  _  ^w 

w  annual  water  discharge  Q 


(23) 


(24) 


87.  To  determine  the  annual  water  discharge  from  DeGray  Lake,  a 
10-year  flow  record  (1966-1976)  for  Glenwood,  Arkansas,  was  used  for  a 
first  approximation.  From  a  previous  study,  which  predicted  the  annual 
inflow  to  DeGray  Lake  (Perrier  and  Ford  1978),  a  correction  of  30  per¬ 
cent  on  an  area  basis  was  required  to  incorporate  the  drainage  area  be¬ 
tween  Glenwood,  Arkansas,  and  the  Highway  84  gage.  The  latter  stream 
gage  was  assumed  to  be  the  inflow  to  DeGray  Lake.  Therefore,  the  annual 
discharge  rate  from  DeGray  Lake  was  estimated  as  590  x  10^  cu  m/yr. 

Using  the  TP  average  annual  loading  in  Table  17,  the  incoming  phosphorus 
concentration  at  steady-state  was 

P  =  annual  loading  .  . 

o°  yearly  water  discharge  "  7 

Table  18  lists  the  estimated  inflow  phosphorus  concentrations  to  DeGray 
Lake  using  the  Highway  84  gage  data.  Equation  21  was  used  to  compute 
the  phosphorus  retention  coefficient  as  follows: 


R  =  0.854  -  0.142  In  q 
exp  s 


where 


q 


s 


Z_  _  Z  X  Q 

T  V 

w  w 


11.04  cu  m/yr  (Highway  84  gage) 


In  q  =  2.4 
s 
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a  ad 


R 

exp 


0.51 


88.  According  to  Figure  1  in  the  Larsen  and  Mercier  (1976)  study, 
the  MRI  and  NES  methods  result  in  DeGray  Lake  being  classified  as  either 
oligotrophic ,  mesotrophic,  or  eutrophic  based  on  the  MRI  range  and  the 
NES  mean  annual  TP  loading.  The  three  instream  methods  place  DeGray 
Lake  in  the  mesotrophic  to  eutrophic  state,  definitely  not  oligotrophic 
(Table  18).  DeGray  Lake  is  most  closely  classified  as  mesotrophic  based 
on  EPA's  interpretation  of  data  obtained  from  the  NES  sampling  program. 
This  application  emphasizes  the  extreme  caution  that  must  be  used  by  CE 
Districts  when  applying  Vollenweider-type  analyses  to  estimate  the  tro¬ 
phic  status  of  a  reservoir  based  on  annual  phosphorus  loading  estimates. 


Table  18 

Estimated  Mean  Annual  Inflow  Phosphorus 
Concentration  to  DeGray  Lake 


Loading  Technique 

Glenwood,  Ark. 

Highway  84 

Trophic  State 

MRI 

22 

29 

Oligotrophic,  mesotrophic 

Range 

11 

- 

275 

11 

- 

275 

or  eutrophic 

NES 

26 

+ 

13* 

34 

+ 

17 

Oligotrophic,  mesotrophic 

or  eutrophic 

Q*C 

42 

+ 

8* 

55 

+ 

11 

Mesotrophic 

QC 

75 

+ 

18* 

97 

+ 

23 

Eutrophic 

FI  @  5  cu  m/sec 

42 

+ 

3* 

55 

+ 

4 

Mesotrophic 

*  ±  1  standard  error,  pg/£  . 
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PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


89.  Conclusions  and  recommendations  based  on  the  results  of  this 
study  are  as  follows: 

a.  The  MRI  procedure  for  suspended  solids  (USLE)  is  based  on 
theoretical  input  requirements  that  are  virtually  impos¬ 
sible  to  quantify  for  any  watershed  perhaps  larger  than  a 
farm  field.  Similarly,  C  (NT)  and  C  (PT)  are  extremely 
important  in  the  loading  calculations  for  nitrogen  and 
phosphorus.  However,  they  are  not  easily  quantified  for 
specific  watersheds  because  of  the  heterogeneous  soil  com¬ 
position  and  land  use  patterns.  This  technique  does  allow 
the  planner  to  become  familiar  with  watershed  physiography 
and  hydrology,  which  enables  him  to  better  interpret  loading 
estimates  for  a  specific  reservoir.  Consequently,  the  MRI 
procedure  is  recommended  for  those  projects  with  little  or 
no  existing  strearoflow  and  water  quality  data,  and  where 

a  projection  of  future  loadings  resulting  from  changing 
land  use  is  required. 

b.  The  NES  procedure  uses  regional  regression  equations  that 
were  developed  using  results  from  a  nationwide  field  study 
that  incorporated  only  1  year  of  data  collection  at 
monthly  intervals  for  each  sampling  location.  The  equa¬ 
tions  are  at  best  applicable  to  geographical  areas  where  a 
significant  number  of  NES  sampling  locations  exist.  There¬ 
fore,  the  NES  method  only  has  limited  application  to  CE 
projects  within  the  vicinity  of  the  watersheds  included  in 
their  survey.  It  is  not  recommended  for  general  use. 

c.  The  average  flow,  average  concentration  Q'-'C  method  is 
only  applicable  to  water  quality  constituents  that  main¬ 
tain  fairly  constant  concentrations  through  the  range  of 
streamflow.  Consequently,  it  is  not  recommended  for  SS , 

TKN ,  and  TP. 

d.  The  applicability  of  the  flow-weighted  concentration  QC 
loading  technique  depends  heavily  on  the  degree  to  which 
the  data  for  a  specific  watershed  are  representative  of 
the  normal  range  of  flow  and  water  quality  concentrations. 
This  technique  is  not  recommended  for  general  use. 

e.  The  flow  interval  (FI)  technique,  which  incorporates 
the  probability  of  occurrence  for  given  flow  and  water 
quality  concentration,  is  a  statistically  valid  method  for 
analyzing  existing  instream  data.  However,  those  water 
quality  constituents  whose  loading  is  dependent  on  season 
and  flow  require  sufficient  data  to  compute  the  short¬ 
term  or  seasonal  loadings  and  sum  them  to  obtain  an  esti¬ 
mate  of  the  annual  loading.  This  technique  is  recommended 
for  analysis  of  instream  data. 
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f.  Field  sampling  programs  to  obtain  an  instream  data  set  for 
estimating  water  quality  loadings  should  include  simulta¬ 
neous  measurement  of  flow  and  quality  over  a  wide  range  of 
streamflow  conditions.  In  addition  to  periodic  baseflow 
sampling,  selected  storm  hydrographs  should  be  sampled 
during  different  seasons  of  the  year.  Results  of  the  Lake 
Erie  wastewater  management  study  (U.  S.  Army  Engineer  Dis¬ 
trict,  Buffalo  1975)  and  previous  investigations  on  the 
Caddo  River  (Westerdahl,  Perrier,  and  Nix  1977,  Westerdahl 
et  al.  1975)  clearly  illustrate  the  importance  of  ele¬ 
vated  flow  sampling.  The  required  sampling  frequency  will 
depend  on  the  hydrologic  characteristics  of  the  watershed 
and  the  magnitude  and  duration  of  the  storm.  Samples  gen¬ 
erally  should  be  taken  at  1-  to  2-hr  intervals  throughout 
the  ascending  limb  and  peak  of  the  hydrograph.  Sampling 
frequency  normally  can  be  reduced  for  the  descending  limb 
of  the  hydrograph.  After  the  hydrograph  is  recorded,  cost 
of  sample  analysis  can  be  reduced  by  selecting  samples 
representative  of  various  portions  of  the  storm  hydrograph 
and  discarding  excess  samples.  Normally,  about  10  samples 
per  storm  are  required  with  emphasis  placed  on  the  ascend¬ 
ing  limb  and  peak  flow  portions  of  the  hydrograph.  Auto¬ 
matic  samplers  that  are  activated  by  changes  in  flow  ve¬ 
locity  or  river  stage  are  recommended  when  extensive  sam¬ 
pling  is  required. 
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